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TECHNICAL MEMORANDUM 


ORIGINAL PAGE m 
OF POOR QUALITY 


AN ANALYTICAL APPROACH TC) THERMAL MODELING OF 
BRIDGMAN-TYPE CRYSTAL GROWTH: 
ONE-DIMENSIONAL ANALYSIS 

Computer Program Users Manual 

INTRODUCTION 


This report is a users guide for a computer program developed to simulate one- 
dimensional heat flow in a rod inserted into a Bridgeman-type directional solidification 
furnace. The thermal model was developed by R. J. Naumann of NASA's Space 
Sciences Laboratory at Marshall Space Flight Center and is included in Appendix B. 
Using this model, a small scientific computer has been applied to the study of several 
situations of interest in the Bridgman-type crystal growth. The computer program 
listing is in Appendix A. 


THERMAL MODEL 


The complete thermal model for the one-dimensional Bridgman-type crystal growth 
is included in Appendix B. A brief description follows, but the complete report 
should be referenced for applications and limitations before attempting to utilize the 
computer program. 

A) The model is a one -dimensional analytical description of heat flow in a trans- 
lating or motionless rod (sample) inside a directional solidification ftirnace with an 
adiabatic zone separating the hot and cold zone. 

B) An optional booster (tickler) heater is included in the hot zone just in 
front of the adiabatic zone. 


HOT ZONE 


BOOSTER 

HEATER ADIABATIC 
ZONE ZONE 


COLD ZONE 


C) The jnodel is applicable to systems whose Biot numbers in the hot and cold 
zones are less than unity. The Biot number for the adiabatic zone is assumed equal 
to zero. 

D) Different thermophysical properties of the sample in the solid and liquid 
state can be accommodated. 

E) The space between the sample and furnace wall can be either vacuum or 
filled with a stagnant or moving fluid. Different heat transfer coefficients in the hot 
and cold zone can be accommodated. 

F) The enclosure of the sample in an ampoule is handled by calculating effec- 
tive conductivities, specific heat, and density as area weighted averages. These in 
turn are used to calculate effective Peclet numbers, heat transfer coefficients, and 
Biot numbers. 
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Reference Tables A- 2 and A- 3 for term definitions. 

G) Three main cases aie treated: 

i) A translating infinite rod in a three-zone (hot, adiabatic, and cold) 
furnace. 

ii) A motionless rod with finite length inserted into the hot zone and 
infinite length in the cold zone. An adiabatic zone is included. 

iii) A motionless infinite rod in a four- zone [hot, booster (optional), 
adiabatic, and cold] furnace. 


COMPUTER PROGRAM 


A flow diagram and complete program listing are given in Appendix A, Key 
program symbols are listed in Table A-1. The program is coded in BASIC language 
for a Hewlett Packard (H-P) 9835 computer with print er/plotter output and occupies 
core storage of approximately 18 K. 

Data are read in via the interactive mode on the keyboard with the required 
units specified. The general input is listed in Table A- 2. 


SELCAS contains an interactive case selection allowing for any one of the 
following. 

Case I: Calculates and plots temperature profiles in a translating infirdte rod. Refer- 
ence pages 5 to 8 of the Thermal Model (Appendix B>. Additional required 
input for tliis case is listed in Table A- 3. The output consists of: 

A) Input 

B) Effective Conductivities 

C) Peclet Numbers 

D) Heat transfer coefficients and Biot numbers 

E) Coefficients used in temperature calculations 

F) Interface position and gradients in the liquid and solid at the interface 

G) Tabulated temperature profiles in the sample. 

NOTE: After the profiles are printed, the temperature at both ends of 

the adiabatic zone (x = ±1/2 zone length) should be checked as to the 
accuracy of tlic input values of Tin and Tsa. An interactive option is 
given to i) correct these input temperatures to agree with the printed 
profiles or ii) continue the program and plot the profiles. These tem- 
peratures (Tla and Tsa) at the outer edges of the adiabatic zone are 
used to calculate the radiative heat transfer coefficients. Page 14 of 
the Thermal Model in Appendix B should be referenced as to their 
importance. 

H) Plot (»f temperature versus distance in the sample. 

NOTH: In output G and H. x will be in the range Lim x him where 

Him is an input value. Fifty one values are printed /plotted. 

t\ise II; Calculates and plots the temperature profiles in a motionless rod with a 

finite length inserted into the lx)t zone. Tliis case can be used to investi- 
gate the end effects which arc not included in Case 1 or 111. Reference 
pages 8 to 10 of the Thermal Model (Appendix B). Additional required input 
is given in Talde A- 4. The output is the same as for Case 1 except no 
Pei’lot numbers !U*e printetl. 

Casj III: Calculates and plots temperature profiles in an infinite motionless rod where 
an optional bixjster (tickler) heater has been auded to the hot end. This 
case can l>e used to investigate interface {wsition and gradients in the steady 
stat(' nKuic with or witlxnit the booster heater. Reference pages 11 and 12 
of the rhormal Model (Appendix B). Additional required input is listed in 
I'ablc A r>. Output has the same format as Case 11. 

Case IV : Calculates t h.e optimum ct>ltl end temperature (Tc) which will center the 

r,, isotherm within 10 cm of the center of the adiabatic zone. The hot 
Melt 

end temperatur<' (I'h) is input. The interface jxjsition Xo (equation (39)1 
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is set equal 0 and equations (38) and (39) solved simultaneously for Tc. 
Reference page 12 in Thermal Model (Appendix B). Temperature profiles are 
then calculated from Case III equations. This case is of interest to crystal 
growers desiring the melt interface in the middle of the adiabatic zone. In 
many cases this position is the optimum for the flattest possible interface. 
Additional required input is the same as Case III (Table A- 5). 

NOTE: The solution for Tc involves an iteration between that equation 
and the Biot number calculation. If the solution cannot be found within 
10 iterat'ons an error message is printed and an option to select another 
case is given. The output has the same format as Case III with the cal- 
culated Tc being printed also. 

Case V: Calculates the optimum length for the adiabatic zone in order to achieve a 

designated liquid gradient at the solid-liquid interface. The booster (tickler) 
heater may be on or off. Equation (38) of the Thermal Model is solved for 
the adiabatic zone length La. Reference page 12 of the Thermal Model 
(Appendix B). Additional required input is listed in Table A 6. The output 
has the same format as Case I parts A, B, and D plus the calculated adiabatic 
zone length. 

NOTE; If the calculated adiabatic zone length is negative this implies the 
designated liquid gradient cannot be achieved. An interactive option is 
given to calculate Case 111 temperature profiles. 

The subroutine BIOTCAL calculates the radiative and conductive heat transfer 
coefficients. These plus the input convective heat transfer coefficients are used to 
compute effective heat transfer coefficients for the bot and cold zones of the furnace. 


(eff) = rad + Hj^ conv + cond 


(eff) = rad + conv + cond 


These effective values are in turn used to calculate the Biot numbers. Reference 
pages 12 to 15 and page 24 of the Thermal Model (Appendix B). If the Biot numbers 
are available as input they should be set in this section. 

TXPIX)T controls the plotting of temperature versus distance in the sample, 

AXE and LABLER subroutines are called from here. These are H P 9835 utility sub- 
routines which draw and label the axes and would probably require changing if the 
program is run on another computer. The "Plotter is" statement should be checked 
for correct set-up even if running on the H-P 9835. 


SAMPLE PROBLEMS 


Each of the five cases was run for a 0. 5 cm radius Mn-Bi/Bi eutectic sample. 
The sample was enclosed in a 0.6-cm radius fused silica ampoide. Furnace configura- 
tion was modeled with the Solidification Experiment System (SES) furnace developed 


by TRW for use on the Space Shuttle. Forced convection with helium gas is used for 
heat extraction in the cold end: stagnant helium gas is in the furnace cavity in the 
l»t end. 

In Cases 111 and IV the booster heater was not turned on, and Tb = Th was 
used to model this. Case V did utilize the bqoster heater. The general input for 
each case was the same and is listed in Table A- 7. The specific case input is listed 
in Table A -8. All output from the five cases follow the Concluding Remarks section. 


CONCLUDING REMARKS 


Current work is being devoted to: 

A) Attempting to evaluate furnace performance by using this model along with 
temperature data from a well characterized and instrumented sample^ to solve for 
effective heat transfer coefficients. 

H) Compiling a users guide to accompany the computer program for R. J. 
Nauniann's "An Analjdical Approach to Thermal Modeling of Bridgman-Type Crystal 
Growth: Two Dimensional Anulysi.s." 

Doth tasks will be the subject of a further report. 


1. Data from NASA's SPAR 9 flight. Mn-Hi/Bi Solidification Experiment. 
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1*0 HDDEL or CKYbTHL OROMTH 
S^hPLE CASE I 


ORIGINAL PAGE IS 
OF POOR QUALITY 


TEMPERATURE OF HOY ZONE tt) 
TEMPERATURE OF lri''ERFACE (C) 
TEMPERATURE OF COl.D ZONE iC> 
CONDUCTIVITY OF L'OUIO (M/c*-Cj 


CONDUCTIVITY OF SO..IO <W/c»-Ci 

CONVECTIVE HEAT TRANSFER COEFF.-HOT END (N/cm 2-C> 
CONVECTIVE HEAT TRANSFER tUEFF.-UOLO tNDfM/cnZ-Cj 
CONDUCTIVITY 0^ GAS - HOT END iM/ca-C) 

CONDUCTIVITY OF GAS - COLD END (M/cm-C) 

CONDUCTIVITY OF AMPOULE (M/cm-C> 

SAMPLE RAOlU' (cm) 

AMPOULE OUTER RADIUS (cm> 

INSIDE RADIUS OF FURNACE MUFFLE IN HOT ZONE (cm) 

INSlOE RADIUS OF FURNACE MUFFLE IN COLO Z0 NE(cmF 

EFFECTIVE EMISSIVITY OF AMPOULE /SAMPLE -HOT END- 
EFFECTIVE EMISSIVITY OF AMPOULE /SAMPLE -COLD ENO- 
EMISSIVITY OF FURNACE -HUT END 
EMISSIVITY OF FURNACE -COLD END 
LENGTH OF ADIABATIC ZONE (cn) 

TEMP. OF SAMPLE AT HOT END OF ADIABATIC ZONE (C) 
TEMP. OF SAMPLE AT COLO END OF AOIABAFIC ZUNE LCJ 
CALCULATED EFFECTIVE CONDUCTIVITY OF LIQUID (H/cw-C) 
CALCULATED EFFECTIVE CONDUCTIVITY OF SOLID (M/cm-C) 


CASE I 


AVERAGE SAMPLE DENSITY 19 / 011 * 3 ) 

AVERAGE AMPOULE DENSITY ( 9 /CM 3 ) 

AVERAGE SAMPIE SPECIFIC HEAT (W-aec/ 9 -C) 
AVERAGE AMPOULE SPECIFIC HEAT (ri-<iec/ 9 -C> 
AMPOULE VELOCITY (Cw/aec/ 

HEAT OF FUSION fW-aec/ 9 ) 

CALCULATED PECLEC HO. FOR THE LIQUID 
CALCULATED PFCLEC HO. FOR THE SUL ID 


CALCULATED 

CALCULATED 

CALCULATED 

CALCULATED 

CALCULATED 

CALCULATED 

CALCULATED 

CAL'^ULATED 

CALCULATED 

CALCULATED 


EMISSIVITY FUNCriOH lFeh)-HOT 
EMISSIVITY FUNCTION (FecZ-CULO 
RADIATIVE MEAT TRANS. COEFF.-HOl (H/cm 2-C) 
RADIATIVE HEAT TRANS. COEFF.-COLD (M/Cm2-0) 
CONDUCTIVE HEAT TRAHS. COEfF. -HUT LM/cm2-C) 
CONDUCTIVE HEAT TRANS. COEF F . -C0L0 (M/cm 2 -C ) 
EFFECTIVE HEAI TRANS. CUE F F . -HO T IN/ c«2-C ) 
EFFECTIVE HEAT TRAHS. COEF F . -C0 LD(M/cm2-C> 
BIOT HO. -HOT EHO- 
BIOT NO. -COLD END- 


4S0 

271.5 

40 

. 124 
.072 
0 

.OJl? 

.003 

0 

.02 

.5 


2.6 

1.77 

. 7695 

.9486 

.3 

.3 

2.03 

327 

133 

9.22222222223E-02 

5.61111111112E-02 


10.05 

2.3 

. 1463 

.493 

.0028 

50.16 

2.491195B4337E-02 
4. 0944406930 7E-02 


.54406/658336 
.541962660261 
4. 11201340730E-03 
4. 72619886225E-04 
3.40965719206E-03 
0 

7.52187059936E-03 
3. 16726198862E-02 
4.893/4713693E-02 
.338677519574 


THE FOLLUHIHG COEFFICIENTS ARE USED IN CALCULATIN<> TEMPERATURE 

A- 1923.32846936 B« 664 . 08424353 ASTAR- 1923. 328469/3 BS T AR« - 1672 . 0 1967207 C- 
91 . 3925254013 D- 122.083895174 


INTERFACE POSITION Xo (cn) 
GRADIENT IN LIQUID AT Xo 
GRADIENT IN SOLID AT Xo 


-. 1/8038124/95 


(de 9 C/cmt 
( de 9 C/c 6 ) 


6 


-69.583803614 
-112. /21895074 


1 


TEMPERATURE PROFILE 


X(c»> 

TtC) 

-5 

435.951856441 

-4.B 

434.341590689 

-4.6 

432.546748534 

-4.4 

430.546172943 

-4.2 

428.31628176 

-4 

425.83078973 

-3.8 

423.06U39e«53 

-3.6 

419.972452034 

-3.4 

416.530550126 

-3.2 

412.694120874 

-3 

408.417941658 

-2.8 

403.651606221 

-2.6 

398.338930504 

-2.4 

392.41 7290353 

-2.2 

385.81688334 

-2 

378.459905943 

-1.8 

370.259636433 

-1.6 

361.1194126 

-1.4 

350. 931 '.92371 

-1.2 

339.575783/46 

-1 

326.92217423 

-.8 

313.6104447 

-.6 

300. 1877146 

. 4 

286.65305831 

- . 2 

273. 00554233 

0 

251.30879766 

. 2 

228.33241068 

. 4 

205.04028899 

.6 

181.4280941 7 

.8 

157.4914277 

1 

133.225830/9 

1.2 

111.352633473 

1.4 

94.6000530217 

1.6 

81 . 7807394768 

1.8 

71.9712178/67 

2 

64.4648320074 

2.2 

58. 7208384589 

2.4 

54.3254526522 

2.6 

50.9620407315 

2.8 

48.3883099483 

3 

46.418854437 

3.2 

44.911 798^42 

3.4 

43. 7585706526 

3.6 

42.8761180033 

3.8 

42.2008465258 

4 

41.684119158 

4.2 

41.2887120047 

4 . 4 

40.9861408103 

4.6 

40. 7546090159 

4 . 8 

40.577437584 

5 

40.4418634768 


ORIGINAL PAGE 13 
OF POOR QUALITY 


Input temperatures at hot and cold 
edge of adiabatic zone were checked 
and program continued. 
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1-D nODLL OF CRVSfAL OROMTH 
SAMPLE CASE I I 


ORIGINAL PaGc i3 
OF POOR QUALITY 


TEMPERATURE OF HOT ZOME (C) 

TEMPERATURE OF INTERFACE CO 
TEMPERATURE OF COLD ZONE CO 
CQHOUCTIVITV OF LIQUID CM/cm-O 
CONDUCT IV! TV OF SOLID CU/cn-O 

CONVECTIVE HEAT TRANSFER COEFF.-HOT END (M/c»2-0 
CONVECTIVE HEAT TRANSFER COEFF.-COLO EH0CM/eii»2-O 
CONDUCTIVITY OF GAS - HOT END (W/e«-0 
CONDUCTIVITY OF GAS - COLO END (M/em-C> 

CONDUCTIVITY OF AMPOULE CM/cm-O 
SAMPLE RADIUS Ccm) 

AMPOULE CUTER RADIUS Ccm) 

INSIDE RADIUS OF FURNACE MUFFLE IN HOT ZONECcm> 

INSIDE RADIUS OF FURNACE MUFFLE IN COLD ZONECcm) 

EFFECTIVE EMISSIVITY OF AMPOULE/SAMPLE-HOT END- 
EFFECTIVE EMISSIVITY OF AMPOULE /SAMPLE -COLD ENO- 
EMISSIVITY OF FURNACE -HOT END 
EMISSIVITY OF FURNACE -COLD END 
LENGTH OF ADIABATIC ZONE CcM> 

TEMP. OF SAMPLE AT HOT END OF ADIABATIC ZONE CO 
TEMP. OF SAMPLE AT COLD END OF AO I ABAT I C ZONE CO 
CALCULATED EFFECTIVE CONDUCTIVITY Of LIQUID CW/c»-0 
CALCULATED EFFECTIVE CONDUCTIVITY OF SOLID CM/cn-O 


450 

271.5 

40 

.124 

.072 

0 

.0312 

.003 

0 

. OZ 

.5 

.G 

2.G 

1.7/ 

.7695 

.9486 

.3 

.3 

2.03 

327 

133 

9.22222222223E-02 

5.G1111111112E-92 


CASE I 1 


BIOT HO. AT END OF ROD IN HOT ZONE .0507 

LENGTH OF AMPOULE INSERTED INTO MOT ZONE (cm) 2.5 


CALCULATED 

CALCULATED 

CALCULATED 

CALCULATED 

CALCULATED 

CALCULATED 

CALCULATED 

CALCULATED 

calculated 

CALCULATED 


EMISSIVITY FUNCTION (Feh)-HOT 
EMISSIVITY FUNCTION (Fec)-COLO 
RADIATIVE HEAT TRANS. COEFF.-HOT (M/cm2-C) 
RADIATIVE HFAT TRANS. COEFF.-COLO (M/Cm2~C> 
CONDUCTIVE MEAT TRANS. COEFF.-HOT iU/emZ-O 
CONDUCTIVE HEAT TRANS. COEFF . -COLDLM/cm2-C) 
EFFECTIVE MEAT TRANS. COEFF . -M0T(M/cm2-C) 
EFFECTIVE HFAT TRANS. COEFF . -C0L0(M/cm2-C) 
GiUT NO. -HOT ENO- 
BIOT NO. -COLO EHO- 


.544967658338 
.541962660261 
4.11201340730E-03 
4.72619886225E-04 
3. 409857192 06E -03 
0 

7.52187059936E-03 

3.16726198862E-02 

4.89374713693C-02 

.338677519574 


THE FOLLOWING COEFFICIENTS ARE DEED IN CALCULATING TEMPERATURE 
B-- 2. 43982292246 C- 78.2141450111 D- 1"-?. 221897134 


INTERFACE POSITION Xo (cm) -.41376185103 

gradient in LIQUID AT Xo (deg C/cm) 

ORAC ENT IN SOLID AT Xo (deg C/cm> 


-65.2762906338 

-107.235/97755 


y 




original page is 

T£|i^ERATUtC 

RftOr lU 

OF POOR QUALITY 

XCcm) 

T<C> 


-2 

359. 134587772 


- 1.92 

356.244814424 


- 1.84 

353. 191883499 


- 1.76 

349.970482121 


- 1.66 

346.575004245 


- 1.6 

342.999540885 


- 1.52 

339.23786964 


- 1.44 

335.283444853 


- 1.36 

331.12938423 


- 1.28 

.' j 26 . 768458865 


- 1.2 

322. 1930 79853 


- 1.12 

317.39528428 


- 1.04 

312 . 36672-385 


-.96 

307.156404538 

Input tefoperatures at hot and cold edge 

-.88 

301.934300647 

of adiabatic zone were checked. They 

-.8 

296. 712196756 

varied from the profile. When the 

-.72 

291.490092865 

Interactive mode allowed for a correction 

-.64 

266.267988975 

the following was made: 

-.56 

281.045885084 

Tla - 310* 

-.48 

275.823781193 

Tsa - 113* 

- . 4 

270.023548834 

The following run uses this new Input. 

-.32 

261 . 440685014 


-.24 

252.857821193 


-.16 

244.274957373 


-.08 

235.692093552 


0 

227. 109229732 


.08 

218.526365912 


. 16 

209.943502091 


.24 

2 01 . 3666.^2 71 


.32 

192 . 7/777445 


.4 

184.19491063 


.48 

1 / 5 . 612046 U 1 


.56 

167.029182989 


.64 

158.446319169 


. 72 

149.863455349 


.8 

141.280591528 


.88 

132 . 697727/08 


.96 

124.114863888 



1.04 

115 . 5 / 746/354 

1.12 

107. 722789311 

1.2 

100 . 66443679 b 

1.28 

94.3775756408 

1..36 

86 . 7261776013 

1 . 44 

83.662122965 

1.52 

79. 1243694444 

1.6 

75. 0582193543 

1.68 

71 .414659502 

1 . 76 

66 . 1497/05755 

1.84 

65. 2241977475 

1.02 

62.6026762923 

2 

60 . 2536065042 


j U MuUll uF L«frsi**i. oyuuiH 
SA«KI tASl II 


TCHPCffATURC OF HOT ZOHf (O 4b0 

TCnPfRATURF OF INTERFACF tc, 2?1.S 

TEWCifFETUHl UF COLO ^OHf (t <0 
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OF POOR QUALITY 

CONOUCTIVITV OF LIQUID (W/c»-C} 

COHDUCTIVITV OF SOLID CM/c«-C) 

CONVECiIVE HEAT TRANSFER COEFF.'HOT END (H/caZ-Cl 
CONVECTIVE HEAT TRANSFER COEFF. -COLD END(M/cm 2-C> 
CONDUCTIVITY OF CAS - HOT END <M/c»-C) 

CONDUCTIVITY OF GAS - COLO END <W/c«-CJ 
CONDUCTIVITY OF AMPOULE (M/cb-C) 

SAMPi E RADIUS (c«) 

AMPOULE OUTER RADIUS (cn) 

INSIDE RADIUS OF FURNACE MUFFLE IH HOT ZONE (cm) 

INSIDE RADIUS OF FURNACE MUFFLE IN COLO ZONE(cm) 

EFFECTIVE EMISSIVITY OF AMPOULE/ SAMPLE-HOT END- 
EFFECTIVE EMISSIVITY OF AMPOULE/SAMPLE-COLD ENO- 
EMISSIVITY OF FURNACE -HOT END 
EMISSIVITY OF FURNACE -COLO END 
LENGTH OF ADIABATIC ZONE (cm) 

TEMP. OF SAMPLE AT HOT END OF ADIABATiC ZONE (C) 

TEMP. OF SAMPLE AT COLD END OF ADIABATIC ZONE (C) 

CALCULATED EFFECTIVE CONDUCTIVITY OF LIQUID (M/cm-C) 
CALCULATED EFFECTIVE CONDUCTIVITY OF SOLID (M/cm-C) 


.124 

.072 

0 

.0312 

.003 

0 

.02 
.5 
.6 
2.G 
1.77 
. 7695 
.9466 
.3 
.3 

2.03 

310 

119 

9.22222222223E-02 

S.61111111112E-02 


CASE II 


AIOT NO. AT END OF ROD IN HOT ZONE .0507 

LENGTH OF AMPOULE INSERTED INTO HOT ZONc (cm) 2.5 


CALCULATED EMISSIVITY FUNuTlON .Feh)-HOT 
CALCULATED EMISSIVITY FUNCTION (Fec)-COLD 
CALCULATED RADIATIVE HEAT TRANS. COEFF. -HOT (M/cm2-C) 
CALCULATED RADIATIVE HEAT TRANS. COEFF. -COLD CM/Cm2-C) 
CALCULATED CONDUCTIVE HEAT TRANS. COEFF. -HOT (H/cm2-C) 
CALCULATED CONDUCTIVE HEAT TRANS. COEFF . -C0LDCM/cm2-C) 
CALCULATED EFFECTIVE HEAT TRANS. COEFF . -H0TCM/cm2-C) 
CALCULATED EFFECTIVE HEAT TRANS. COEFF . -C0L0(M/cm2-C) 
CALCULATED BIOT NO. -HOT END- 
CALCULATED BIOT NO. -COLD END- 


.544067658336 
.541962660261 
4. 042388287S2E-03 
4.56648445440E-04 
3. 40985719206E-03 
0 

7.45224S47958E-03 
3. 165664844S4E-02 
4 . 84844886623C-02 
.338506735851 


THE FOLLOWING CuEFFlCIENTS ARE USED IN CALCULATING TEMPERATURE 
B--2. 4844750218C C- 78.0336904002 0- 132.577142304 


INTERFACE POSITION Xo (cm) -. 419? 13415946 

GRADIENT IN LIQUID AT Xo (deg C7c») 

GRADIENT IN SOLID AT X. Cdeg C/cm7 


-65. 109271623 
-107.011278113 


TEMPERATURE 

PROFILE 

X(cn>) 

TCC) 

-2 

358.604688279 

-1.92 

355. 71687064 

-1.84 

352 . 666495634 

-1.76 

349.448303983 

-1.68 

346. 056747072 

-1.6 

342.485977374 

1.52 

338. 72983839 

-1.44 

334 . >81854 

■1.36 

330.635217328 

-1.28 

326. 28277899 

-1.2 

321. 717034768 

-1.12 

316.930112676 

-1.04 

311 . 91375939 
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-.96 

306.716731518 

-.88 

301 . 5079897 ^ 

-.8 

296.299248058 

-.72 

291.090506329 

-.64 

285.881764599 

-.56 

280.673022869 

-.48 

275.464281139 

-.4 

269.45464893 

-.32 

260.893746681 

-.24 

252.332844432 

-.16 

243.771942183 

-.08 

235.211039934 

0 

226.650137685 

.08 

218.089235436 

.16 

209.528333187 

.24 

290.967430938 

.32 

192.406528689 

.4 

183.84562644 

.48 

175.284724191 

.56 

166. 723821942 

.64 

158.162919692 

.72 

149.602017444 

.8 

141.041115195 

.88 

132.480212945 

.96 

123.919310697 

1.04 

115. 403748068 

1.12 

107.56899419 

1.2 

100.54830286 

>-* 

• 

ro 

CD 

94.2570897066 

1.36 

88.6195590011 

1.44 

83.5677905018 

1.52 

79. 0409211469 

1.6 

74,984411797 

1.68 

71.3493901474 

1.76 

68.0920619269 

1.84 

65. 1731832614 

1.92 

62.5575878739 

2 

60.2137634087 


Input Tia and Tsa were checked and 
program allowed to continue. 
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SRMPLE CASE II 

TEMP(C) 



3 
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1-0 MODEL OF CRYSTAL OROMTH 
SAMPLE CAbt 111 


TEMPERATURE OF HOT 20HE ID 

temperature OF IMTERFACE <C) 

TEMPERATURE OF COLO ZOME <C) 
conductivity of liquid CM/etn-C) 

CONDUCTIVITY OF SOLID (M/cm-C) 

CONVECTIVE HEAT TRANSFER COEFF.-HOT ENO (H/e«>2-C> 
CONVECTIVE HEAT TRANSFER COEFF.-COLO LH0 (H/cm 2-C> 
CONDUCTIVITY OF GAS - HOT EHD 
CONDUCTIVITY OF GAS - COLO EHD (H/ca-D 
CONDUCTIVITY OF AMPOULE (M/ca-D 
SAMPLE RADIUS (ea> 

AMPOULE OUTER RADIUS Cca) 

INSIDE RADIUS OF FURNACE MUFFLE IN HOT TONECca) 
INSIDE RADIUS OF FURNACE MUFFLE IN COLO ZONECca) 
EFFECTIVE EMISSIVITY OF AMPOULE/SAMPLE-HOT EHD- 
EFFECTIVE EMISSIVITY OF AMPOULE/SAMPLE-COLD EHO- 
EMISSIVITY OF FURNACE -HOT ENO 
EMISSIVITY OF FURNACE -COLD END 
length OF ADIABATIC ZONE (ca) 

TEMP. OF SAMPLE AT HOT END OF ADIABATIC ZONE (C) 
TEMP. OF SAMPLE AT COLD FHO OF ADIABATIC ZONE <C> 
CALCULATED EFFECTIVE CONDUCTIVITY OF LIQUID IW/ca-C) 
CALCULATED EFFECTIVE CONDUCTIVITY OF SOLID tM/ca-C) 


450 
271. S 

4« 

.124 

. 0/2 

0 

.0312 

.003 

0 

.02 
.5 
.6 
2.6 
1.77 
. 7695 
.9486 
.3 
.3 

2.03 

327 

133 

9.22222222223E-02 

5.61111111112E-0Z 


CASE I I I 

TEMPERATURE OF TICKLER HEATER (CJ 450 

LENGTH OF TICKLER HEATER icm) 0 


CALCULATED EMISSIVITY FUNCTION (FehJ-MOT 
CALCULATED EMISSIVITY FUNCTION (Fec7-COLD 
CALCULATED RADIATIVE HEAT TRAHS. COEFF.-HOT iM/ca2-C> 
CALCULATED RADIATIVE HEAT TRAHS. COEFF.-COLO (M/Ca2-C> 
CALCULATED CONDUCTIVE HEAT TRAHS. COEFF.-HOT CM/ca2-C> 
CALCULATEO CONDUCTIVE HEAT TRANS. CDEFF . -C0LDlM/ca2-C) 
CALCULATED EFFECTIVE HEAT TRAHS. CuEFI . -H0TtM/ca2-C7 
CALCULATEO EFFECTIVE HEAT TRAHS. COEFF . -C0LD(W/ca2-C> 
CALCULATED BIOT NO. -HOT EHU- 
CALCULATEO BIOT HO. -COLD EHO- 


.544067658338 
.541962660261 
4. 11201340730E-03 
4. 7261988822SE-04 
3. 409857192 06E-03 
0 

7.521870S9936E 03 
3. 16726198862E-02 
4.89.K/4713693E-02 
.338677519574 


THE FOLLOWING COEFFICIENTS ARE USED IN CALCULATING TEMPERATURE 
A--130. 957228179 B- 0 C« 81.8170404875 0- 130.95/2281/9 


INTERFACE POSITION Xo (ca) -.318741494102 

GRADIENT IN I IQUIO AT Xo (des C/ca) 

GRADIENT IN SOLID AT Xo (dC 9 C/ca7 


-68.283218687 

-112.227864377 


TEMPERATURE PROFILE 


X (ca) 
5 

-4.8 

4.S 

-4.4 

4.2 

-4 


T(C) 

433.804166315 
431.802021864 
429.801795101 
<27.581713854 
425. U 761 2864 
422. 3826 7C 759 


*3.8 

419 

.347118092 

-3.6 

415 

.977913321 

-3.A 

412 

.238383152 

-3.2 

408 

. 087823344 

-3 

403 

.48105567 

-2.8 

398 

.367936141 

-2.6 

392 

.692809199 

-2.< 

386 

.393901934 

-2.2 

379 

.40265166 

-2 

371 

.642959655 

-1.8 

363 

. 030362802 

-1.6 

353 

.47111426/ 

-l.< 

342 

. 861 163034 

-1.2 

331 

.085021341 

-1 

318 

. 018523541 

-.8 

304 

.361879804 

-.6 

290 

. 705236066 

-.4 

277 

. 048592329 

-.2 

258 

.173895704 

0 

235 

. 728322829 

.2 

213 

.282749954 

.4 

190 

.837177078 

.6 

168 

.391604203 

.8 

145 

.946031328 

1 

123 

.500458452 

i.2 

103 

. 479837122 

1.4 

88. 

2494159627 

1.6 

76. 

6731586969 

1.8 

67. 

8743388284 

2 

61 . 

1865787593 

2.2 

56. 

1033817634 

2.4 

52. 

2397725065 

2.6 

49. 

3031409932 

2.8 

47. 

0710817782 

3 

45. 

3745501169 

3.2 

44. 

0850593821 

3.4 

43. 

10495014J 

3.6 

42. 

3599939413 

3.0 

41. 

7937716055 

4 

41 . 

363400353 

4.2 

41 . 

03628606/2 

4.4 

40. 

7876547858 

4.6 

40. 

5986764478 

4.8 

40. 

4550388008 

S 

40. 

3458634643 


OR*G*?<Al- i 1'^ 
OF POOR QUkUTY 


Input Tla and Tsa checked. The 
following corrections were made: 

Tla = 318® 

Tsa = 123® 

Next case uses these values. 


10 MUOtL OF CkVSlHL OSrUWTM 
SAMPLE CAbE 111 


TEMPfRATURf Of HOT ZOHE (C) <50 

TEMPERATURE OF INTERFACE iC) 2/1.5 

TEMPERATURE OF COLD ZONE (C) <0 

CUNOUCTIVITV ul LIQUID (W/c«C) . 12< 

CONDUCT I VI TV OF SOLID iW/CBi-L) .0/2 

CONVECTIVE HEAT TRANSFER COEfF. HOT END (W/c«2-CJ 0 

CONVECTIVE HEAT TRANSFER COEFF. COLD ENDtW/c«2-C> .0312 

CONDUCTIVITY OF GAS - HOT END (W/c»-C> .003 

CONDUCTIVITY OF GAS COID LNO iU/ct»-CT 0 

CONDUCTIVITY OF AMPOULE <W/c*-C> .02 

SAMPLE RADIUS <c«) .5 
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Af«»OULE OUTER RADIUS (cn) 

INSIDE RADIUS OF FURNACE MUFFLE IN HOT ZONE(cm) 
INSIDE RADIUS OF FURNACE FHIFFLE IN COLD ZONE(c«> 
EFFECTIVE EMISSIVITV OF AMPOULE/SAMRLE-NOT tNO- 
EFFECTIVE EMISSIVITV OF AHROULE/SAMRLE-COLD END* 
EMISSIVITV OF FURNACE -HOT ENO 
EMISSIVITV OF FURNACE -COLO ENO 
LENGTH OF ADIABATIC ZONE (cmJ 

TEMP. OF SAMPLE AT HOT END OF ADIAB.'ITIC ZONE CO 
TEMP. OF SAMPLE AT COLD END OF ADIABATIC ZONE fO 
CALCULATED EFFECTIVE CONOUCTIVirV OF LIQUID CM/cm-O 
CALCULATED EFFECTIVE COHDUCTIVITV OF SOLID CH/cw-O 


.G 
2.6 
1.77 
. 7695 
.9486 
.3 
.3 

2.03 

318 

123 

9.22222222223C-02 

5.61111111112E-02 


CASE III 

TEMPERATURE OF TICKLER HEATER CO 450 

LENGTH OF TICKLER HEATER (cfn) 0 


CALCULATED EMISSIVITV FUNCTION (Feh^-HOT 
CALCULATED EMISSIVITV FUNCTION (FeOCOLO 
CALCULATED RADIATIVE HEAT TRANS. COEFF.-HOT (M/c»2-C> 
CALCULATED RADIATIVE HEAT TRANS. COtFF.-COLO lM/C«2-0 
CALCULATED CONDUCTIVE HEAt TRANS. CULFF.-HOT (.N/CM2-0 
CALCULATED CONDUCTIVE HEAT TRANS. COEFF . -C0L0(W/c«2-O 
CALCULATED EFFECTIVE HEAT TRANS. COEFF . -HUUW/ci»2-0 
CALCULATED EFFECTIVE HEAT TRANS. COEFF . -COLD tH/c(il2-0 
CALCULATED BIOT NO. -HOT END- 
CALCULATED BIOT NO. -COLD ENO- 


.544067658338 
.541962660261 
4.U7496086894E-03 
4.61152580853E-04 
3. 409657192 06E-03 
0 

.007484618061 
3. 16611525809E-02 
4 . 86964 0666 1 9E - 0 2 
. 338554896884 


THE FOLLOWING COEFFICIENTS ARE USED IN CALCULATING TEMPERATURE 
A- -131. 14432183 8 - 0 C- 81.7466779081 0- 131.14432183 


INTERFACE POSITION Xo (cm) -.320/58831802 

GRADIENT IN LIQUID AT Xo (809 C/cn) 

GRADIENT IN SOLID AT Xo Cde 9 C/cm> 


-68.2121434732 
-112. 111047689 


TEMPERATURE PROFILE 


X(CM) 

T(C) 

-5 

433. 496393694 

-4.8 

431.687113483 

-4.6 

429.679482996 

-4.4 

427.451757213 

-4.2 

424.97980/222 

-4 

422. 236858874 

-3.8 

419. 193202786 

-3.6 

415.815872551 

-3. 4 

412. 068287681 

-3.2 

407.909857386 

-3 

403. 295540935 

2.8 

398.175359809 

-2.6 

392. 493856374 

■2.4 

386. 18949321 

2.2 

379. 193986592 

2 

371 .431566887 

-1.3 

362.818157877 

- .6 

353.260466137 

-1.4 

342.654970534 

-1.2 

330.886800966 

1 

31 7. 832496018 


-.8 304. 190067323 
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. b 

290.547638629 

.4 

276.905209934 

.2 

257.961600849 

0 

235.539391311 

.2 

213.117181773 

.4 

190.694972235 

.6 

168.272762698 

.8 

145.85055316 

1 

123.428343622 

TTT" 

103.428158508 

1.4 

88.2125309016 

1.6 

76.6469434178 

1.8 

67.8557967573 

2 

61.1735370149 

2.2 

56.0942684076 

2 4 

52.2334532674 

2.6 

49.2987997379 

2.8 

47.0681331486 

3 

45.3725757751 

3.2 

44.0837615609 

3.4 

43. 1041178725 

3.6 

42.3594785403 

3.8 

41.7934689374 

4 

41.3632380097 

4.2 

41.0362141393 

4.4 

40.7876392345 

4.6 

40.5986943626 

4.8 

40.4550750193 

5 

40.3459081731 


ORIGINAL PAGE EGT 
OF POOR QUALITY 


Input T1a and Tsa were checked 
and program allowed to continue. 



ORIGINAL P/vj: i 
OF POOR QUAUi V 


l-D nOOEL or CKVbTAL oROMTH 
bAHPLE CASE IV 


TEMPERATURE OF HOT 20HE (C) 
temperature of interface cC) 

TEMPERATURE OF COLO ZONE IC> 

CONDUCTIVITY OF LlOUIO CM/cm-C) 

CONDUCTIVITY Or SOLID (W/c»-C) 

CONVECTIVE HEAT IRANSFER COEFf -HOI END 
CONVECTIVE NEAT TRANSFER COEFF. -COLO EN0 LM/cm 2-C) 
CONDUCTIVITY OF GAS - MOT END 
CONDUCTIVITY OF GAS - COLO ENO (M/cm-C) 

CONDUCTIVITY OF AMPOULE (M/cm-C) 

SAMPLE RADIUS (ca) 

AMPOULE OUTER RADIUS (c«> 

INSIDE RADIUS OF FURNACE MUFFLE IN HUT ZUNElcia) 

INSIDE RADIUS OF FURNACE MUFFLE IN lOLO ^UNELcn^ 

EFFECTIVE EMISSIVITY OF AMPOULE / SAMPLE-HOT ENO- 
EFFECTIVE EMISSIVITY OF AMPOULE/ SAMPLE-COLO ENO- 
EMISSIVITV OF FURNACE -HOI ENO 
EMISSIVITY OF FURNACE -COLO END 
LENGTH OF ADIABATIC ZONE (cm) 

TEMP. OF SAMPLE AT HOT END OF ADIABATIC ZONE (C) 
TEMP. OF SAMPLE AT COLD END OF ADIABAIIC ZONE (C) 
CAICULATED EFFECTIVE CONDUCTIVITY OF LIQUID (H/cm-C) 
CALCULATED EFFECTIVE CONDUCTIVITY OF SUClD (M/cm-C> 


4S0 
271.5 
40 
. 124 
. 0/2 
0 

.0312 

.003 

0 

.02 
.5 
.6 
2.6 
1.77 
. 7695 
.9486 
.3 
.3 

2.03 

327 

133 

9.22222222223E-02 

5.61111111112E-02 


CASE IV Tc CALCULAIIOH 


CALCULATED Tc(C) 

TEMPERATURE OF TICELER HEATER IC) 
LENGTH OF IICELER HEATER (cm) 


97.270796601 

4S0 

1.8 


CALCULATED 
CALCULATED 
CALCULATED 
CALCULATED 
CAl CULATED 
CALCULATED 
CALCULATED 
CAL CULATED 
CALCULATED 
CALCULATED 


EMISSIVITY FUNCTIUH (Feh)-HOT 
EMISSIVITY FUNCTION (Fec)-COLO 
RADIATIVE HEAT TRAMS. (.UEFF.-HOT lH/c»2-C) 
RADIATIVE HEAT TRANS. COEFF. -COLO i.M/Ci»2-C) 
CONDUCTIVE HEAT TRANS. CUEfF.-MOT (M/cm 2-C) 
CONDUCTIVE HEAT TRAHS. COtf F . -COLO(W/c»2-C) 
EFFECTIVE HEAT TRANS. COEF F . -H0T(H7ciii2-C) 
EFFECTIVE HEAT TRAHS. COEF f . - C0 LD(W/cm2 -C) 
BIOT HO. -HOT LHD- 
BIOT NO. COLD EHO- 


.544067658336 
.541962660261 
4. 11201340730C-03 
6. 71110029012E-04 
3.40985719206C-03 
0 

7.52187059936E-03 
.031671110029 
4.89374713693E -02 
.340799986428 


THE FQLIOWING COEFFICllNTS ARE USED IN CALCULATING lEMPERATUKE 
A- 73.0062345341 B« 0 C- 72 . 697 7544 731 0* 45.6621749258 


INTERTACE POSITION Xo (cm) 

GRADIENT IN LIQUID AT Xo (deq C/CM) 
GRADIENT IN SOLID AT Xo (dog C/c«) 


-2. 65152250032E-08 


-60.8622235947 

-100.030981353 


TEMPERATURf PROF 1 1 1 
X(cm) T(Ci 


4 . A 
4 . 6 

4 . 4 


T35. 36605949 
433 . ) 79771434 
431 . 996926581 
43!). 01812494b 
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-4,2 

427.F»21821b4 7 

-4 

426.384111442 

-3.8 

422.678462024 

-3.6 

419.675421188 

-3.4 

416.34230178 

-3.2 

412.642823287 

-3 

<08.636717428 

-2.8 

< 03.9/9289848 

-2.6 

IM. 920933621 

-2.4 

393 . 3065892 78 

-2.2 

387.07514551 

-2 

300.158773977 

-1.8 

372.40?!<?5963 

-1.6 

363.961838068 

-1.4 

354.504972427 

-1.2 

344.008657485 

-1 

332.362221869 

-.8 

320. 18977714 

-.6 

308.017332421 

-.4 

296.844887702 

-.2 

283.672442983 

0 

271.499997148 

. 2 

261.4938008/7 

.4 

231.487604607 

.6 

211.481408336 

.8 

191.4/6212065 

1 

171.469015794 

1.2 

153.630440521 

1.4 

140.071582928 

1.6 

129. 774680279 

1.8 

121.954979681 

2 

116. 016522809 

2.2 

111.506/24406 

2.4 

108. 081881668 

2.6 

105.480978732 

2.8 

103.50579431 / 

3 

102.U05/94(S67 

J.2 

100.86666106 

3.4 

100.001577274 

3.6 

99.3446131116 

3.8 

98.8456996908 

4 

98.4E>66135 724 

4.2 

98. 1 790/8898 

4.4 

97.9605666884 

4.6 

97. 794623582 

4.8 

97.6686026396 

S 

97.5/28994901 


Input TIa and Tsa were checked, 
following corrections were made 

Tla * 332° 

Tsa - 171° 

Next run uses this Input. 


l-U MODJ L OF COYSIAl OKOMTH 
iAHPlt CAbt IV 


TEMPERA rURE OF HOT ZOHt (C) A50 

TEMPERATURE OF INTERFACE ICT 

TEMPERATURE OF COLO ZONE (C) 47 . Z FO 7lM>601 

CONDUCTIVITY OF LIQUID (W/c»-C) .124 

CONDUCTIVITY OF SOLID (W/cm-t} .072 

CONVECTIVE MEAT TRANSFER COEFF . MOT END <M/c»2-C> 0 

CONVECTIVE HEAT TRANSFER COEFF. -COLO END(M/ca2-C7 .0312 

CONDUCTIVITY OF GAS - HOT END tW/caCy .003 

CONDUCTIVITY OF GAS - COLD END (W/ca-C) 0 
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CONDUCTIVITY Of NNPOUCE iH/ca-C) 

SAMPLE RADIUS (CM) 

AI^OULE OUTER RADIUS (cm) 

INSIDE RADIUS OF FURNACE MUFFLE IN HOT ^ONEicm) 

INSIDE RADIUS OF FURNACE MUFFLE IN tOLU ZONliCMi 

EFFECTIVE EMISSIVITV Of AMPOULE /SAMPLE -NOT ENO- 
EFFECTIVE EMISSIVITV OF AMPOULE /SAMPLE -COLO END- 
EMISSIVITV OF FURNACE -MOT ENU 
EMISSIVITV OF FURNACE -COLD END 
LENGTH OF AOIARATIC ZONE (c») 

TEFtf*. OF SAMPLE AT HOT END OF AOIABAliC ZONE IC) 
TEMP. OF SAMPLE AT CULO END OF ADIABATIC ZONE LC) 
CALCULATED EFFECTIVE CONDUCTIVITY OF LtOUiO (H/cm-C) 
CALCULATED EFFECTIVE CONDUCTIVITY OF SOLID (H/cm-C> 


.02 

.S 

.€> 

2.6 

1.77 

. 7695 

.9486 

.3 

.3 

2.03 

332 

171 

9.22222222223C-02 

5.61111111112E-02 


CASE IV fc CALCULATION 

calculated Tc(C> 

TEMPERATURE OF TICKLER HEATER (C) 

LENGTH OF TICKLER HEATER (cm) 


97.175483728 

450 

1.8 


CALCULATED EMISSIVITV FUNtriOH (Feh)-MOT 
CALCULATED EMISSIVITV FUNCTION (Fec)-COLD 
CALCULATED RADIATIVE HEAT TRANS. COEFF.-HOT (H/cm2-C> 
calculated radiative heat TRANS. COEFF.-COLO (U/Cm2-C> 
CALCULATED CONDUCTIVE HEAT TRANS. COEFF . -HOI LM/cmZ-C) 
calculated CONDUCTIVE HEAT TRANS. COE F F . -CULDim/cm2-C> 
calculated EFFECTIVE HEAT TRANS. CUEF F . -HUT (N/CM2-CJ 
CALCULATED EFFECTIVE HEAT TRANS. COEFF . -COLOTM/cmZ-C) 
calculated BIOT NO. -HOT END- 
CALCULATED BIOT HO. -COLD ENO- 


.544067658338 

.541962660261 

4.13278618099E-03 

7.25332379524C-04 

3.40985719206E-03 

0 

7. 54264 3373056 '03 
3.1 9253323 795E- 02 
4.90/26195355E-02 
.341379/9178 


THE FOLLOWING COEFFICIENTS ARE USEU IN CALCULATING TEMPERAlURE 
A* -72.9241498065 B- 0 C> 72 . 701502031 3 U- 45.5813093083 


INTERFACE POSITION Xo (cm) 

GRADIENT IN LIQUID AT Xo (d «9 C/cm) 
GRADIENT IN SOLID AT Xo Tde 9 C/cm) 


1.282808633546-08 


60.9171141793 

100.121197562 


TEMPFRAIURE PROFILE 


X(cm) 

I vC) 

5 

435.43485571 

4.8 

433.831604654 

4.6 

432. 0518/653 

4 . '» 

430.076245714 

4.2 

427.883148313 

4 

425.448646799 

.T .8 

422. /46168734 

3.6 

419. 746216734 

3.4 

416.416046499 

3.2 

412. 719309425 

3 

408.615655841 

2.8 

404.060294622 

2.6 

399. 003504268 

2 . 4 

393.390090216 

2 . 2 

38 7. lSH/8'2306 

? 

380 . 24 1566429 

I . 8 

372.562941349 

1 . 6 

364. 039095407 


MP ▼ CSJ CM ^ (A CD 


< \m 


-1.4 3S4. 576991338 

-1.2 344.Q73350843 

-1 332.41>11496i 

3503355^24 

308.1)50269289 
295.668846453 
283.683423617 
271.500000781 
251.475761772 
231.451522259 
211.42/282747 
191.403043235 
1 171.378803722 


1.2 

153.525834/88 

1.4 

139.959200417 

1.6 

129.658802454 

1.8 

121.838277827 

2 

115.900584382 

2.2 

111.392420913 

2.4 

107.969620/8^ 

2.6 

105.3708772«- 

2.8 

107,3977944<ii 

3 

101.899741426 

3.2 

100.762352294 

3.4 

99.8987953961 

3.6 

99.2431442249 

3.8 

98. 7453444165 

4 

98.36/3924503 

4.2 

98. 0804343474 

4.4 

97.8625628/11 

4.6 

97.6971450637 

4.8 

97.571552437 

5 

97.4761968756 


owGiNM- 


Input T1« and Tsa chicked and program 
allowed to continue. 
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1 -u nooCL ur uKvsirfL urumth 

crse V 


ORiGiNAL r . : : 
OF POOR QUAik 


TEMRERRTURE OF MOT 20HE (O 
TEnRCRRTURE OF INTERFACE (C> 

TCnfERATURE OF COLO 20NE KC) 

CONOOCTIVITY OF LtOUIO (H/cm-C) 
conductivity of SOLID tN/cm-CJ 

CONVECTIVE HEAT TRANSFER COEFF.-HOT END (M/c*2-C) 
CONVECTIVE HEAT TRANSFER COCFF.-COLO EN0(i4/c«2-C J 
CONDUCiiVITY OF GAS ' HOT END (W/c«-CJ 
CONDUCTIVITY OF GAS - COLD END ^M/e*-C> 

CONOOCTIVITY OF AMPOULE (W/cni-C) 

SAMPLE RADIUS <c«0 
AMPOULE OUTER RADIUS i,cm> 

INSIDE RADIUS OF FURNACE MUFFLE IN HOT 20NElc«> 
INSIDE RADIUS OF FURNACE miFTLE IN COLO ZONELcw) 
EFFECTIVE EMISSIVKY OF AMPOULE/SAMPlE-NOT END- 
EFFECTING EMISSIVITY OF AMPOULE /SAMPLE- COLD END- 
EMISS1VITV OF FURNACE -NOT END 
EMISSIVITY OF FURNACE -COLD END 
LENGTH OF ADIABATIC ZONE (c«) 

TEMP. OF SAMPLE AT HOT END OF ADIABATIC 2UHE (C) 
TEMP. OF SAI«>LE AT COLD EHO OF ADIABATIC ZONE <c; 
CALCULATED EFFECTIVE CONDUCTIVITY OF LIOUID (H/cb-C) 
CALCULATED EFFECTIVE CONDUCTIVITY OF SOLID 


4S0 
2/1. S 
-10 

. 124 
. 0/2 
0 

.0312 

.003 

0 

. 02 

.5 

.6 

2. G 

1.77 

. /GIFS 

.94UG 

.3 

.3 

2.03 

327 

133 

9.22222222223E-02 

S.61111111112E-02 


CASE V calculates UPTIMUM La FuR GIVEN GRADIENT 

>00 

1 .0 


TEMPERATURE OF TICKLER HEATER <C) 
LENGTH OF TICKLER HEATER (ca) 


CALCULATED EMISSIVITY FUNCTION (Fch)-HOT 
CALCULATED EMISSIVITY FUNCTION TFec>-C0L0 
CALCULATED RADIATIVE HEAT TRANS. COEFF.-HOT (H/cb2-CJ 
CALCULATED RADIATIVE HEAT TRANS. COEFF.-CGLD (H/C«2-C) 
calculated CONDUCTIVE HEAT TRAHS. COEFF.-HOI <M/ca2-CI 
CALCULATED CONDUCTIVE HEAT TRANS. COEF F . -COLOLM/ c*2 - C > 

calculated effective heat Trans, coeff. -huT(m/cm2-c; 

CALCULATED EFFECTIVE HEAT TRANS. CUEE F . -COLDTN/ cb2 -C ) 
calculated BIOT NO. -HOT ENO- 
CALCULATED BIOT NO. -COLT/ ENO- 


.S440G76SH33B 
.S419b2bb0261 
4.11201340730E-03 
4. ;261‘MMG22SC-04 
3.4090S71920GE-03 
0 

;.S2l8/0S993GC-03 
3. 16;261Be8G2F-D2 
4.N93/4;i3b93C-02 
.338G77S19S74 


FOR LiaUlD GRADIENT * 95 (C/cai AD. ZONE LENGTH- 1 . 4G832 /’.^1313 lea) 

NOTE I IF AD. ZONE LEHoTH IS NEGA I I VE I NHUI GRADIEHI CAHNUT BE ACNIEVCO WITH 
GIVEN PARAMETERg. 


NOTEi iHl AO. ZONE lENGTH HAS ULEN CHANGED IN KRUG. 

I’ROF llES 


RUN CASE I 1 I FOR TEMP. 


APPENDIX A 


TABLES, FLOW CHARTS, AND COMPUTER PROGRAM LISTINGS 



TABLES 


ORIGINAL PAGS m 
OF POOR QUALITY 


Computer Program Symbol 

Thermal Model Notation 

A 

A 

Af 

af 

Ain 

inner radius of ampoule 

Alfa 

a 

Ao 

a 

o 

Aprime 

A* 

Astar 

a* 

B 

B 

Beta 

B 

Bh , Be 

B. 

i 

Bhstar 

8* 

H 

B prime 

B* 

Bstar 

p :ic 

C 

C 

Cps, Cpa 

speeifie lieat 

D 

D 

Dhf 


Efhot, Efcold 

c, furnaee 

Ehot, Ecold 

c, sample/ampoule 

Feh, Fee 

Fe 

FI 

F(L) 

G1 


Glxo 


Gs 

G 

s 

Gsxo 

G {X ) 
s o 

Heondh, Heondc 

^eond 

Heonvh , Heonvc 

O 

conv 

Hradh, Hrade 

H ^ 

"ad 

Hstarh, Hstare 

H* 

Kgh, Kgc 

^G 





Computer Program Symbol 

Thermal Model Notation 

Kl, Kls 


Ks, Kss 


Ka i 

cxjnductivity of ampoule 

L 

L 

La 


Lb, Lbt 


Lf 

-f 

PI, Ps 

P 

e 

Rhos, Rhoa 

r> 

T 

T 

Tb 

T 

B 

Tc 

T 

Th 


Tla 

Tl(-Xj) 

Tm 

Tm 

To 

To 

Tostar 

To 

Tsa 

T,(Xj) 

U 

u 

X 

X 

Xo 


XI 


X2 

^2 

X3 

-3 







TABLE A- 2. GENERAL INPUT FOR ALL CASES 


Variable 

Units 

Description 

Id$ 


18 Character Run Identification 

Th 

°C 

Temperature of Furnace Muffle - Hot Zone 

Tm 

°C 

Melt Temperature of Sample 

Tc 

°C 

Temperature of Furnace Muffle - Cold Zone 

Kls 

W/cm-°C 

Conductivity of Sample Liquid 

Kss 

W7cm-°C 

Conductivity of Sample Solid 

Hconvh 

wycm^-°c 

Heat Transfer Coefficient due to Convection - Hot Zone 

Hconvc 

W7cm^-°C 

Heat Transfer Coefficient due to Convection - Cold Zone 

Kgh 

W /cm - °C 

Conductivity of Gas in Furnace Cavity - Hot Zone 

Kgc 

W7cm-°C 

Conductivity of Gas in Furnace Cavity - Cold Zone 

Ka 

W7cm-°C 

Conductivity of Ampoule 

Ain 

cm 

Sample Radius 

Ao 

cm 

Ampoule outer Radius 

Afh 

cm 

Inside Radius of Furnace Muffle - Hot Zone 

Afc 

cm 

Inside Radius of Furnace M\iffle - Cold Zone 

Ehot 


Effective Emissivity of Ampoule7 Sample - Hot Zone 

Ecold 


Effective Emissivity of Ampoule/Sample - Cold Zone 

Efhot 


Emissivity of Furnace Muffle - Hot Zone 

Efcold 


Emissivity of Furnace Muffle - Cold Zone 

Lat 

cm 

Total Length of the Adiabatic Zone 

Tla 

°C 

Sample Temperature at Edge of Adiabatic Zone - Hot Side 

Tsa 

°C 

Sample Temperature at Edge of Adiabatic Zone - Cold Side 

Lim 

cm 

-Lim < X < Lim is range to Calculate Temperature 
Profiles. Must be an integer (plotter requirement). 


TABLE A- 3. SPECIFIC INPUT FOR CASE I 


Variable 

Units 

Description 

Rhos 

g/cm^ 

Average Sample Density 

Rhoa 

g/cm^ 

Average Ampoule Density 

Cps 

W-sec/ 

g-°C 

Average Sample Specific Heat 

Cpa 

W-sec/ 

g-°C 

Average Ampoule Specific Heat 

U 

cm /sec 

Ampoule Pull Rate 

Dhf 

W-sec/g 

Heat of Fusion 
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TABLE A- 4. SPECIFIC INPUT FOR CASE II 


Variable 

Units 

Description 

Bhstar 


Biot Number p.t End of Sample - Hot Zone 

Lc 

cm 

Length of the Ampoule Inserted into the Hot Zone 


TABLE i 

V-5. SPECIFIC INPUT' FOR CASE III AND IV 

Variable 

Units 

Description 

Tb 

°C 

Temperature of Booster (Tickler) Heater 

Lbt 

cm 

Length of the Booster (Tickler) Heater 


TABLE A- 6. SPECIFIC INPUT FOR CASE V 


Variable 

Units 

Description 

Tb 

°C 

Temperature of Booster (Tickler) Heater 

Lbt 

cm 

Length of Booster (Tickler) Heater 

Git 

- °C / cm 

Gradient in Liquid at which to Optimize Adiabatic Zone 



Length. Note: Must be negative. 
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TABLE A-7. GENERAL INPUT FOR CASES I-V 
(Reference Table A- 2 for definitions and units) 


ld$ 

Sample Case I (see Note 1} 

Th 

450 

Tm 

271.5 

Tc 

40 

Kls 

0. 124 

Kss 

0.072 

Hconvh 

0 (see Note 2) 

Hoonvc 

0.0312 (see Note 3) 

Kgh 

0.003 

Kgc 

0 (see Note 4) 

Ka 

0. 02 (see Note 5) 

Ain 

0.5 

Ao 

0. 6 (see Note 6) 

Afh 

2.6 

Afc 

1.77 

Ehot 

0.7695 

Ecold 

0.9486 

Efhot 

0.3 

Efcold 

0.3 

Lat 

2.03 

Tla 

327 (see Note 7) 

Tsa 

133 (see Note 8) 

Lim 

5 (see Note 9) 


NOTES 

1. I was replaced by II. Ill, IV. and V for the five cases respectively. 

2. Heat transfer in the hot end was assumed to be dominated by conduction and 
radiation since the space between the sample and furnace muffle contained stagnant 
helium gas. A value of Hconvh = 0 was used here to model that assumption. 

3. The convective heat transfer coefficient in the cold zone was estimated from the 
formula in the Thermal Model (page 15. Appendix B). 

4. The conductive heat transfer of the gas was included in the empirical formula 
for forced convection. Thermal Model (page 15. Appendix B). 

5. A fused silica ampoule was modeled and an average conductivity value used. 

Enter 0 here if no ampoule is included. 
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TABLE A- 7. (Continued) 


NOTES (Continued) 

6. Enter Ao = Ain if no ampoule is used. 

7. If this value is not known, a good first guess is Tla = (Th + Tm)/2. The 
interactive mode will allow for a correction, if necessary after the temperature 
profiles are calculated. 

8. See note 7 and use Tsa = (Tc + Tm)/2 as the first gfuess. 

9. Lim was set = 2 (integer required) in Case II because only 2.5 cm was inserted 
into the hot 7X)ne. 


TABLE A 8. SPECIFIC INPUT FOR CASES I-V 


Case I 


Rhos = 10. 05 
Rhoa = 2. 3 
Cps = 0.1463 
Cpa = 0. 493 
U = 0.0028 
Dhf = 50.16 

Case 11 


Bhstar = 0.0507 (see Note 1) 



Lc = 2.5 

Case III 


T^ = 450 (see Note 2) 
Lbt = 0 (see Note 3) 

Case IV 


Tb = 450 (see Note 2) 
Lbt = 1.8 (see Note 3) 

NOTE: An input for Tc = 40 was entered when general input was 

called for. This case will calculate a new value for Tc and run the 
temperature profiles with this new value. 

Case V 


Tb = 500 



Lbt = 1.8 



Git - -85 (see Note 4) 


NOTES 

1. Case III was run and the printed temperature profile was used to find T at 

X Lc - 2.5. This T was then read in as input for Tla in Case III, and 

the I’alculated Biot number for the hot zone was used as Bhstar input for 
Case 11 . 

2. If the booster heater is not used, this value should be set equal Tl;. 

3. If Iho booster heater is not used (i.e.. Tb = Th), Lbt ruay be any number > 0. 

4. This must be a negative number. 


















































ORIGINAL PA3H IS 
OF POOR QUALITY 


COMPUTER PROGRAM LISTING 


jt 

20 

50 

60 

70 

99 


90 

PIH Tfmi 

109 

Begin > 

llA 

C«ne'*0 

1^0 

Pd«»I »l 

130 

INPUT 

lAO 

INPUT 

ISO 

IMPUT 

'.GO 

INPUT 

170 

INPUT 

IBO 

1 HPU f 

190 

INPUT 

203 

INPUT 

210 

IMPUT 

220 

IHPUT 

2.30 

INPUT 

240 

INPUT 

250 

INPUT 

260 

INPUT 

2 70 

I HPU T 

280 

IHPUT 

290 

INPUT 

50 0 

input 

.310 

1 nput 

320 

INPUT 

330 

IHPUT 

IS 

NOT KHOI 

.340 

IMPUt 

3S0 

input 

3€0 

8r9 1 n2 1 

3T'C 

PRIn ( 

38. 

PRINT 

390 

PRINT 

4C-0 

PPlHl 

4 1 (5 

PRINT 

4.'0 

PWtHT 

■I.IO 

PRINT 

440 

F P I H I 

45 0 

PRINT 

4fc0 

PRIM 

4 7<3 

P'.i'jHI 

■ifO 

PR 1 h: 

490 

PR I H I 

*.0 0 

PPIWT 

* 1 1 

PB 1 N 1 

0 

?'W INI 

',30 

PRIHI 

'40 

PR ! M I 

%-,o 

PRINl 

SfcO 

PRINT 

s n 

P8IKT 

SHO 

PBIHl 

•.98 

PBIHI 


PRINT 

Oi u 

PR 1 H I 

tj *' ft 

FR i N 1 


AH HHAiVIICAl ro tHiHIMk. nVUlLINV 

Of B«IOOIWH-Typt t«V>rHL UfFOMIH 

OHr-biniLNSIOHAl. AHALYSIS 

Hl&fkTIVI « 1'^. IM MOl /0«l. Of fUDHAGC 
»>0 IH Hli^DU; Of ODiABAtU. /OHE 
POSITIVE * IS IH COl.O ^OM£ Of f OHHACt 


Count.«r-0 • I UWAT I OH f OUNICB 
•CASI SfltCIIOH COHfWOl. 

0 !CAS£f COHTBOl VAtlABt.1 

' BUN lOENTIf ICATIOH (IB CHABACtEBST * . I dl 
•TtHf. Of ;iOr ?OHt ( Otfe. l,)*,Th 
‘TfHA. Of Mill IHILBfACt tOtO. C)*,l» 

•TCHP (If COtO ^OHC fOEO C>",lc 
•COHDOCtlvlIV Of LIUUID m/ir«-C)*.M» 

•COHOUCT IVI TV Of SOUP («/«:•- C I ",* »» 

‘COHVEik IVE HEAT TRANSfER COEff, - HOT EHD (M/caZ'C> * .HeaHvh 
•COHVECTIVt HEAT TRAHbflR COEff. -COLO EHO (M/ etn;; -C > * .Hcenwc 
•COMDWCTlVI IV Of GAS -HOT EMO ( M/c«-C)*,K9H 
•C9HOOCTIVI TV Of OAS COLO EHO- IM/ CH- C ) * , K^C 
•COHDUCTlVITV Of AWOULE ( W/c*-C)*,f« 

•SAW'LE RAOlUS ( c»)*.Aif* 

•Af^ROUlE OUUR RADIUS tc*i)*.Ao 

•INSIDE RADIOS Of fURNACE HUffLl IH HOT ^OHt (e«>*,Afh 
“IHStOE RADIOS Of f ORMACE HUf f I E IH THE COLD ^OHE lCM>*,Afc 
•fffECTiVf EHISSlVirV Of AW'UUtE /SArtPLE -HOT EH(l-*,ENol 
•EffECnVE EHlSSIVlIV Of AMPOOLE/SAMPlE-LOtO CHD'.Ecold 
•EHISSIVITV Of lURHACE WAlI HOI CHO’.EfNot 
•EHISSlVirv Uf TURHALE hall COLO EHO'.Efcold 
•ifHOTH Of AOlABATIt 2^UHl ^cm)•.Lmt 

•ft HP. Of SAHPIE A I MOl tHU Of AO I ABA: 1C ^OMt (DEO C)*.fl» ‘If 
WM OUJSS A VAlUf.CHfLf Af TLR PROfUE IS CAL. AHO LORRECT 
•TEHR. Uf sAHKE AT COLD t nO Uf fOIABAtlC /UHE (DEO Cf.Tia 
•St PJ TJ.'IP. PROfllt HILL Bf cal. f OR -LIHdS^LlH. IHPUl L IH(CH> 


PRIMIER IS .’,6 


'L'RlHItR 


I 0 


HODl I 
t Id* 


01 CRYSIAL OROHTH^ 


• rrHPI RflUlRL 

• UHf I RrtlORl 

• UHPLRAUiRf 
•COHU'JCI IVI IV 


Of 

Of 

Uf 

Of 


•COHOOI T t VI f V Of 
'COMVtCTIVE HI A I 
•COHVI Cnvt HI A I 
•COMPUCIIVIIV 01 
•COHOUCT IVI TV Of 
•COHDUtt I Vi TV if 
•SAHfLi RASmiS ‘ 


HOI /UHt (Cl 
IHILRIACt (C) 

COLD /UHl vLI 
LiUUlD (W/CH C I 
SOI ID (H/ca-Cl 
IRAHSIlR culm . 

IRANSf I R CUEM . 

OAS MOl EAD (H/CH-CI 
GAS COLO IHD (H/ca-CI 
AMPOUll 1H/C« C> 

• ) 


HOT IMP (H/caiif 
CULP EHPlH/car 


•AHPOOll (lUtfR RAPIDS (c»i 

•IHMDf RADiUS U( I URHAtl HUM I E IH HOI ZuNMch) 

•IHSIPI RADIUS Uf f ORHALt HUf I t I IH LOLO /DHE(i;hI 

•flllClIVl EHISSIVIIV Of AHROULL /SAMPLE -HOT EHO- 
•IffECMVt IHISSiVlIV Of AHPOUlE/SAHPLE-tUlO EHD ^ 
•EHISSIVIfV Of fURMACt HO! IMP 

•IHiSSiVlIV Of fURNAll LOl D EHD 

*1>H(>IH Uf AOIABAIIL /ONI (CHI 

• HHP Uf SAHPII A1 HU I t HD Uf AD I AI>A I I C /OHC (C / 
•II HP Ul sAHPI I A I (.in D I HU Uf ADIABATIC /'JHL (Cl 


• , Ih 

• . Ih 

• , 1 c 

• .*!■ 

• ,KS* 

• .MeoifvN 

• . H,. onvc 

• .* 9 C 

• ,Rd 

• ,Al« 

• ,Ao 

• , Aih 

• ,Alc 

• ,CNot 
, Eco I d 
• ,t tHol 
* , L I CO id 

• ,L«l 

• .lid 

• , l»d 


TEMP. 

* ,L IH 



‘ * r* - r’ ' 

ixtvji*-’ 4 . .'LJ, ^ 

OF POOR QUALilY 


630 

640 • NORMALIZE BV AO AND CALCULATE EFFECTIVE CONDUCT I VI TIES 

650 La*L«t/Ao 

660 Xl-La/2 'HALF LENGTH OF AD. ZONE 

670 Xinc«ZaLiM/50 >X INCREMENT FOR TEMP. PROFILES-SET FOR Si POINTS 
680 »c>-' -< 'INITIAL X(ca7 AT WHICH TO START TEMP. PROFILES 

690 Xi-i(l9aAin“2/Ao*2*Kaa{l-Aln“2/Ao*2> (EFFECTIVE CONOOCTIVITV-CIOUIO 
700 K9*K99aAkn*2/Ao*2*Ka«(l-Alii*2/Ao*2) 'EFFECTIVE CONOUCTIVITY-SOLiO 

710 PRINT 'CAL'jULATED EFFECTIVE CONOUCTIVITV OF LIOUID (M/ea-C>>,Ki 
720 PRINT 'CALCULMED EFFECTIVE CONDUCTIVITY OF SOLID tH/c«i-C} *,K9 

730 1 ------ 

740 • SELCAS 

750 • folloning Allows various case selections 

760 ' 


770 

SslCB9 1 

PRIH IS 16 'CRT 



780 

PRINT ‘THE 

FOL. t ?IG RvN OPTIONS ARE AVAILABLE* 

790 

PRINT 




800 

PRINT ‘1. 

CASE 1* 



810 

PRINT ‘2. 

CASE II" 



820 

PRINT ‘3. 

CASE III* 



830 

PRINT ‘4. 

CALCULATE tCASE lll> OPTIMUM 

Tc 

TO CENTER THE HClT ISOTHERM* 

840 

PRINT ‘5. 

CALCULATE (CASE III) OPTIMUM 

La 

FOR A GIVEN GRADIENT* 

850 

PRINT 




860 

INPUT ‘INPUT 1,2, 3. 4, OR 5‘,Ca9e 



870 

PRINTER IS 7,6 'PRINTER 



880 

IT Caac-1 

THEN Caael 



890 

IF Caae-2 

THEN Caae2 



900 

IF Caae-3 

THEN Caac3 



910 

IF Caae-4 

THEN Caac4 



920 

IF Caae-5 

THEN Caae5 




930 PRINTER IS 16 'CRT 

940 PRINT ‘ERROR IN LAST INPUT--ROH OPTION HOST -1,2. 3, 4, OP 5* 

950 PRINT ‘PRESS: CONTINUE BAR FUR ANOTHER CHANCE TO INPUT RUN OPTION* 

960 PAUSE 
970 GOTO Selca; 

980 • 

990 ' 

1000 X------- ---------------- ----------- 

1010 • BIOICAL 

1020 ' CALCULATES CONDUCTIVE AND RADIATIVE HEAT TRANSFER COEFF. AND BIOT NOS. 
1030 ' 

1040 Biotcal: PRINTER IS 7,6 'PRINIER 
1050 Feh-1/ Cl /Ehot«Ao/Af had /Ef hot-1 J) 

1060 Fec-1 /(l/Ecold«Ao/AfcaCl/Efcold-l)) 

1070 Hradh-5.67E- 12aFeha(25/ 12a cTh-2/3)*3* 13/ 12a(Th«273)*2acT la‘2 73) ♦ 7/l?aOh*27 
3)a(Tla-273'>‘*2*l/4a(Tla-2/3)*3) 

1080 Hradc-5.67E-12*FecaC25/12aCTc»273)*3‘13/12a».Tc«2/3)*2aLT9a«273)‘7/12aCTc*27 


3)a(T9a» 
1090 Hco 
1100 Hco 


1110 
1120 
1130 
1140 
1150 
1160 
il/D 
1180 PRi 
1190 PRI 


Hat 
Hat 
Bh* 
Be • 
Alf 
Bet 
IF 


1200 

121v. 

1220 

1230 


PRI 

PRI 

PRI 

PRI 


273)“2d/4aCTaa-273)*3) 
ndh-Kgh/Ao/LOOCAfh/AO 
ndc -Kgc / A o/LOGC AT c / Ao) 
arh-Hconvh*Hradh*Hcondh 
arc - Hconve-Hrade -Hcondc 
HatarhaAo/X 1 
HatarcaAo/Ka 
a-SQR(?aBc) 
a*SQRC2aBh) 

Paaal-l THEN RETURN 
NT 
NT 

NT ‘CALCULATED ENISSlVIfV FUNCTION (Feh)-HOT *,Feh 

NT ‘CALCULATED EHISSIVITY FUNCTION (FeO-COLO '.Fee 

NT ‘CALCULATED RADIAMVE HEAT TRANS. COEFF. -HOT CW/cm2-C) ‘ ,Hradh 
NT ‘CALCULATED RAOIAIIVE HEAI TRANS. COEFF. -COLO CW/C»2-C) * ,Hradc 


'EFFECTIVE HEAT TRANSFER CUEFF.-HOT END 
•EFFECTIVE HEAT IRANSVER COEFF. - COLD END 
'BIOT NO. -HOT END 
'BIOT NO. -COLD END 


ALLOWS FOR TC ITERATION WITHOUT PRIHT-CASE4 


4| * 
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1240 PRINT *CALCULATEO CONDUCTIVE HEAT TRANS. COEFF.-HOT (M/e«2-C} * .Heondh 
1250 PRINT 'CALCULATED CONDUCTIVE HEAT TRANS. COEFF . -COLDCM/c*2-C> * .Hcondc 
1260 PRINT 'CALCULATED EFFECTIVE HEAT TRANS. COCFF . -HOT(H/cn 2*C) '.HstarN 
1270 PRINT 'CALCULATED EFFECTIVE HEAT TRANS. COEFF . -C0LDlM/c*2-C> ' ,Hit*rc 
1280 PRINT 'CALCULATED BIOT NO. -HOT END- ',BN 

1290 PRINT 'CALCULATED BIOT NO. -COLO END- '.Be 

1300 RETURN 
1310 ' 

1320 •--- - - 

1330 ' 

1340 • CASE I 

1350 > FOLLOWING CALCULATES TEHP. PROFILES IN A SAMPLE AS IT IS BEING TRANSLATED 
1360 ' AT A VELOCITY 0 

1370 ' 

1380 Caaeli PRINT 
1390 PRINT 

1400 PRINT ' CASE f 

141C PRINT 
1420 PRINT 

1430 INPUT 'AVERAGE SAMPLE DENSITY L 9 /ci> 3 > ' . Rhoa 
1440 INPUT 'AVERAGE Al>^OULE DENSITY (. 9 / c*3> * . Rhoa 
1450 INPUT 'AVERAGE SAMPLE SPECIFIC MEAT (W-aec/ 9 -CJ ' ,Cpa 
1460 INPUT 'AVERAGE AMPOULE SPECIFIC MEAT VW-aec/g-CJ ' ,Cpa 


1470 INPUT 'AMPOULE VELOCITY tci»/acc7',U 
1480 INPUT 'HEAT OF FUSION (W- aec/g) ' . DNf 

1490 PRINT 'AVERAGE SAMPLE DENSITY C 9 /ciii 3 > * ,Rhoa 

1500 PRINT 'AVERAGE AMPOULE DENSITY (.g/ca3> ' >Khoa 

1510 PRINT 'AVERAGE SAMPLE SPECIFIC HEAT iW-aec/g-CJ *.Cpa 

1520 PRINT 'AVERAGE AMPOULE SPECIF lu HEAT (W-aec/g-CJ ' ,Cpa 

1530 PRINT 'AMPOULE VELOCITY (Ca/aec) ',U 

1540 PRINT 'HEAT OF FUSION IW-aec/g) ' ,Uhf 


1550 Rceff “RhoaaCpaaAi n*2/ Ao“2»RhoaaCpaaL 1 -Ai n“2/ Ao“2) 

1560 Dhfef f -0hf»Alt>“2/Ao“2 

1570 Rhoef f »RhoaaA 1 n*2/Ao*2»Rhoaa( 1 -Ai n“2/ Ao“2 ) 

1580 Pl*RceffaAoaU/Kl 
1590 Pa-IClaPI/Ka 

1600 PRINT 'CALCULATED PECLEC NO. FOR THE LIQUID *,PI 

1610 PRINT 'CALCULATED PECLEC NO. FOR THE SOLID * ,Pa 

1620 GOSUB Biotcal 

1630 Lf >RhoeffaAo«UaOhfef f 

1640 Aatar-CSQRC4»Alfa“2»Pa*2)-Pa)/2 

1650 Batar«CSQR(4aBeta“2'Pl“2) •PlJ/2 

1660 Laatl* CPaaKla(Th-T«>a(Xl* 2 / 2 *Xl/Aatar 7 'PlBK 3 aLTa-Tc)a(Xl* 2 / 2 *Xl/ 9 atar)*Lfa(X 
l-l/Aatar)acxial/Batar))/(Kia(Th- I»> ♦Kaat T«- Tc ) J 

1670 LaM2*LfaO/Batar-l/Aatar}*PaaKaa(Ta-Tc)a«.xi'l/Batar)-PiaKlH(Th-Ta)«(Xl'l/Aa 

tar ) 

1680 Lan»2-lai»2/ (XlaCTh-Ta) ♦XaBCTa-TcJ ) 

1690 La»3*(Lf *PaaCK taiTh- Tm) / 2rKaaC Tia- Tc ) ) *Pia (XaaCTia- Tc J /2*Ma< Th-Tm) > > / <K latTh 
-T»)*Xaa(T»-Tc)) 

1700 OxO>Laal-Laa2axO-LaM3aXO'‘2 

1710 XOp’XO'DxO 

1720 GOSUB Reaidue 

1730 ' - --------- 

1740 Raal-Rea 
1750 X0p-X0p*.01 
1760 GOSUB Reaidue 
1770 Rea2*Rea 

1780 XO-XOp*. 0iaRea2/(Real-Rea2> 

1790 IF ABStRealX . 00001 THEN Pop • TEST FOR CONVERGENCE 

1800 XOp-XO 

1810 GOSUB Reaidue 

1820 GOTO 1740 

1830 ' ---------- 


40 


Ffiv-'i kji 

Or POOR QUALITY 


1840 Residue I '.rea-KlvPlBt fm-ThlMEXPlP l>X0p) / lEXPl -P1«X1 J»CP1 /Bstar- 1 J *eXP(Pl«X0p 
)) 'Ll 

1850 Rres'RsaPaatTc- T»>«ExPlPsaX0p)/ (CMPCPs»Xt)«(P9/A«tar*ll -eXP(PsHXOp) > 

1860 Rea*Rres-Lres 

1870 RETURN 

1880 ■ -------- 

1890 Pop: • XO HAS CONVERGED 

1900 B-(T«-Th)/(EXP(-Pl»Xl)dlPl/Bst«r-l)*£XP(Pl»XO)) 

1910 Bpr lute- (Tc -Ta)/(E XP(P9«X1)M(P9/A9 tar *1) -EXP (PsAXOn 

1920 A-B«EXP(-PiaXl)«(Pi/B9tar-l)»Th 

1930 Apr lme>-Bprl»eaEXP(P9AXl)a(P9/A9tar *l)*Tc 

1940 r- -PsABpr tae«EXP(P9*Xl 1 / Astar 

1950 0--PlAB«EXPt-Pl»Xl)/39tar 

1960 Gl'PlABAEXPCPlaxO) 

1970 09*P9«Bpr iae>EXP(P9«X0> 

1980 PRINT 

1990 PRINT 'THE FOLLOHINC COEFFICIENIS ARE USED IN CALCULATING TEHPERATORE* 

2000 PRINT *A-‘ lAi *8-* i8i ‘ASTAR** ,ApriMei‘BSTAR«* tBpr laei ‘C*‘ |Ct *D-‘ |0 
2010 FOR I-l TO 51 ' CALCULATE 51 PT. TEW. PROFILE 

2020 X-Xc/Ao 'NORHALIZE 

2030 IF X<-X1 then T-Th-DAEXP(B9taraix*Xl)> 

2040 IF CX>-X1> AND (X<X0) THEN T-A»BAEXP(P1»X) 

2050 IF (X>X0) AND CX<X1> THEN T«Apr iaie*Bpr imeaEXPLPsax) 

2060 IF X>X1 THEN T *Tc »CaEXP< - Aataratx -XI > > 

2070 Xarray(!)-Xc 
2080 Tarray(l)-T 
2090 Xc*Xc*Xinc 
2100 NEXT I 

2110 G05UB Prt2 'PRINT OUTPUT 

2120 GOTO Begin >PROGRAH SHOULD NURMALLV NEVER EXECUTE THIS 
2130 ' 

2140 ' ---- 

2150 I CASE II 

2160 ' INVESTIGATION OF END EFFECTS-SAWLE IS INSERTED A FINITE DISTANCE LC IN 
2170 ' THE HOT ZONE. SAWLE IS NGilONLtSS. 

2180 ' 

2190 CaacZi PRINT 

2200 PRINT • CASE I I * 

2210 PRINT 
2220 PRINT 

2230 INPUT 'BIOT NO. (B»T AT END OF ROD IN THE HOT ZONE ■ .flhstar 

2240 INPUT ‘LENGTH OF THE AWOULE WHICH IS INSERTED :NTO HOT ZONE Cca)‘,Lc 

2250 PRINT ‘BIOT NO. AT END OF RUD IN HOT ZONE * ,Bh9tar 

2260 PRINT ‘LENGTH OF AWOULE INSERTED INTO HOT ZONE Urn) ‘ ,Lc 

2270 GOSUB Biotcal 

2280 L-Lc/Ao 'NORMALIZE 

2290 FI>((l>Bh9tar/Beta)»EXP(BetaAL)*vl-8h9tar/6eta)«EXP(-BetaAL))/C(l*Bh9tar/Be 
ta)aEXP(BetaaL) (1 -6h9tar/Beta)MEXP(-Bela«L> ) 

230 0 GI • - (Th Tm»K9/X|B(T*- Tc ) ) / ( La ‘ 1 / A 1 f 9 *F 1/Beta) 

2310 Gs-tCllGl/Xs 

2320 X0-(Kia(Th-r»)a(XlM/Alfe)-X9a(Ta-Tc)«tXl‘H/Beta))/LXia(Th-T»)*X9aCTia-TcJ) 
2330 TO-Tm-GlaXO 
2340 TOstar -Tie-GsaXO 
2350 C- Gs/Alfa 

2360 B’GlaEXPC -BetaBL) /Belt'/ ( ( 1 *Bh9 tar /Be ta J / C 1 -Bha tar /Beta)aEXP(Be taMLL *La> ) -E < 
P( BetalLL La) )) 

237 0 D>-Gl/Beta/(1 - ( 1 - Bha tar /Bet a ) / < 1 *Bh9 tar /Beta)aEXP( -2aBe taWL ) > 

2380 PRINT 

2390 PRINT ‘THE FOLLOWING COEFFICIENTS ARE USED IN CALCULATING TEWERATURE* 

2400 PRINT *B- ‘ jBi ‘C-‘ iCi ‘D- ‘ iO 

2410 FOR I-l TO 51 'CAL. 51 PT . lEW. PROFILE 
2420 X-Xc.'Ao 'NORMALIZE ^ 

2430 IF X< XI THEN T • Th DaEXP (Be taai X -XI ) > *BaEXP( -BetaBl X -X 1 ) ) 
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2440 IF AND <¥<X0> THEN T-T0*OlBX POOR QUALITY 

2450 IF CX>X0> AND (X<X1) THEN T-TCatar *0«aX 
2460 IF X>X1 THEN T>Tc*CaEXP(-Alf«»CX-Xl>> 

2470 X«rr«y(I)-Xc 
248ii^ TarrayCD'T 
2490 Xc-Xc»Xinc 
2500 NEXT I 

2510 GOSUB Prt2 'PRINT OUTPUT 

2520 GOTO Begin 'PROG. HILL NOT NORHALLV EXECUTE THIS STEP 
2530 • 

2540 ■ 

2550 - 

2560 > CnaE III 

2570 • STEADY STATE CALCULATION HlTH AH OPTIONAL TICKLER HEATER 
2580 • 

2590 Cas«3i IF Caae*4 THEN GOTO 2700 'INPUT HAS ALREADY BEEN SUBMITTED 

2600 X1>Lb/2 ■ 1/2 OF AO. ZONE LENGTH 

2610 PRINT 

2£20 PRINT • CASE 111* 

2630 PRINT 

2640 INPUT 'TEMP. OF TICKLER HEATER lUEG. C>*,Tb 
2650 INPUT *LEHGTH OF T.CKLER HEATER lea) *,Lb. 

2660 PRINT 'TEMPERATURE OF TICKLER HEATER <C) '.Tb 

2670 PRINT 'LENGTH OF TICKLER HEATER Cca) * ,Lbt 

2680 GOSUB Blotcal 

2690 • CALCULATES GRADIENT AND XO POSITION FRim CENTER LINE 

2700 Lb'Lbt/Ao 'NORMALIZE LENGTH TO AMPOULE OUTER RADIUS 

2710 Lb2-Lb/2 

2720 X2-Xl»Lb2 

2730 X3-X2*Lb2 

2740 Ealfa>EXP(AlfaaLb2) 

2750 Ebeta*EXP(8etaaLb2) 

2760 B--CTb-Th)/2/Ebela 

2770 Gl*-CTb-Tm'Ka/Kia(Ta-Tc)»(Th-Tb>BEXPl-BetaHLb>>/(La*l/Bela*l/Alfa) 

2780 X0--(Tb-Ta«<Th- Tb)aExP(-Betaa2aLb2>-Ks/Kiat r»- Tc>*Giaa /Beta 1/ Alfa) )/2/Gl 
2790 IF Caae<>4 THEN GOTO Caae3p 'SKIP FULLOMING TEST IF NO^ IN CASE4 
2800 IF A8S(X0)>lE-6 THEN GOTO Caae4p 'ALLOHS FOR TC ITERATION 
2810 > 

2820 • XO IS WITHIN TOLERANCE FOR CENTER OF AO. ZONE ON CASE 4 
2830 PRINT 'CALCULATED Tc(C) '.Tc 

2840 PRINT 'TEMPERATURE OF TICKLER HEATER (C) • , lb 

2850 PRINT 'LENGTH OF TICKLER HEATER lea) ' ,Lbt 

2860 Paasl-0 'RESET CASE TO ALLOW FOR PRINT IN BIOTCAL 
2870 GOSUB Biotca) 

2880 ' 

2890 Caap3p> Ga-KlaGI/Ka 
2900 C*-Ga/Alfa 

2910 A*(B/Ebeta*Gl/Beta)/Ebeta 
2920 0 -A/Ebeta*BBEbeta 
2930 T0*Ta-GiaX0 
2940 TOa-Ta&BaxO 
2950 PRINT 

2960 PRINT 'IHE FOLLOWING COEFFICIENTS ARE USED IN CALCULATING TEMPERATURE* 

2970 PRINT 'A-',A»'B*',Bi'C-*,Ci'D-'iO 
2980 FOR I -1 TO 51 

2990 X*Xc/Ao 'CAL. 51 PT . TEMP. PROFILE 
3000 IF XC-X3 THEN T Th-DaEXP(Betaa(.X*X3) ) 

3010 IF (X>-X3) AND (X<--X1) THEN T • lb' AaEXPlBetaa ( X«X2) ) *BaEXP( -8etaa(X*X2) ) 
3020 IF (X'-Xl) AMD (X<-X0) THEN T-TO'GiaX 
3030 IF (X>X0) AND (X^-Xt) THEN l-TOa'GaaX 
3040 IF X>X1 THEN T-Tc‘CatxP(-Altaa(X-XD) 

3050 Xarrayl I ) -Xc 
3060 TarrayC I ) *T 





3070 Xc>Xc*Xlnc C; < *- ^ • 

3080 NEXT I 
3090 COSUB Prt2 

3100 GOTO Begin 'PPOGRAH HILL HOT HORHBLLV EXECUTE THIS STEP 
3il0 I-------------- - 


3120 • CASE IV 

3130 • CALCULATES Tc HHICH MILL CENTER THE Tneii ISOTHERN 
3140 • TEMPERATURE CALCULATIONS USES CASE 111 EQ. 

3150 ' 

3160 C.«9e4t PRINT 
3170 PPIHT 

3180 PRINT • CASE IV Tc CALCULATION* 

3190 PRINT 

3200 Pasal'l I CONTROL FOR FIRST PASS IN 8I0TCAL 
3210 INPUT *TEMP. OF TICKLER HEATER (Oeg C)*,Tb 
3220 INPUT "LEHCTH OF TICKLER HEATER tcaJ'.Lbt 
3230 GOSUB Blotcal 
3240 Lb'Lbt/Ao 'NORMALIZE 

3250 Ca9e4pi • ENTER CASE4 FOR ITERATION LOOP 
3260 Counter-Counter*! 'ITERATION COUNTER 
3270 IF Counter<10 THEN GOTO Case4pp 

3280 PRINT * ^ : TO CENTER Ta ISOTHERM IN MIDDLE OF AO. ZG.**: MAS NOT FOUND AFTE 

R 10 TRYS* 

3290 PRINT *Tc-* ,Tc, *X0-* ,X0 
3300 GOTO Selcas 
3310 • 

3320 Ca9e4ppi 

3330 Tc-Tma<La*l/8eta*l/Alfa)/(La/2*l/Beta}-KI/K9a(La/2*l/Aifa7«iTb-rM*CTh-Tb)«iE 


CASE V 

CALCULATES OPIIMUM LENGTH OF ADIABATIg ZONE (LA> TO ACHIEVE 
DESIGNATED LIQUID GRADIENT 


PRINT 


CASE V 


CALCULATES OPTIMUM La i OR GIVEN GRADIENT* 


XP(-BetB«LbT*K9/KlaTa)/<La/2*l/Beta> 

3340 GOSUB Biotcal 'RECALCULATE ALFA AND BETA FOR NEXT ITERATION 
3350 GOTO Ca9e3 
3360 
3370 
3380 
3390 
3400 
3410 

3420 CaneSi 
3430 PRINT 
3440 PRINT * 

3450 PRINT 

3460 INPUT *TEMP. OF TICKLER HEATER lOeg C)*,Tb 
3470 INPUT *LENGTH OF TICKLER HEATER (ca)*,Lbt 

3480 INPUT *GRADIEHT IK LIQUID CMIHUS C/Ca) AT MHICH TO OPTIMIZE AO. ZONE LENGTH 
* ,GIt 

3490 PRINT TEMPERATURE OF TICKLER HEATER (C) * , Tb 

3500 PRINT 'LENCTH OF TICKLER HEATER Cca) *,Lbt 

3510 Lb*Lbt/Ao 'NORMALIZE 
3520 Gl-Glt»Ao 'NORMALIZE 
3530 GOSUB Biotcal 

3540 La-( -Tb*Ta* vTb-Th)«EXPC-BetaaLb)-K9/Kia(Ta-Tc7-Giau/Alfa*l/8eta>)/01 
3550 PRINT 
3560 PRINT 

3570 PRINT "FOR LIQUID GRADIENT - * ,GI t , * (C/ca) AD. ZONE LENGTH- *» La«Aoi * (cm) * 
3580 PRINT 

3590 PRINT *N(JTE: IF AD. ZONE LENGTH IS HEGAT I VE- - I HPUT GRADIENT CANNOT BE AC 

HIEVED MITH* 

3600 PRINT * 

3610 PRINT 
3620 PRINT 
3630 PRINT "NOTE! 

I f OR TEMP. * 


GIVEN PARAMETERS.* 


THE AO. ZONE length HAS BEEN CHANGED IN PROG. RUH CASE II 


3640 PRINT • 


profiles* 
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f'RTZ 

OUTPUT AND PROORAH CONTROL 


Prt2i ' OUTPUT AND PROORAH CO 

PRINT 

PRINT 

PRINT 'INTERFACE POSITION Xo (c«> 
PRINT 'GRADIENT IN LIQUID AT Xo 
PRINT 'GRADIENT IN SOLID AT Xo 
PRINT 
PRINT 

PRINT ' TEHPERATURE PROFILE* 
PRINT ' X(c«) TC 


Cde9 C/cm) 
ideg C/cm) 


‘ ,X0*Ae 

* .Gl/Ao 
* .C*/Ao 


3760 PRINT ' TEHPERATURE PROFILE* 

3770 PRINT * X(c«) TCO* 

3780 FOR I-l TO 51 ‘PRINT 51 PT. fEHP. PROFILE 
3790 PRINT Xarray(I) .TarrayU) 

3800 NEXT I 

3810 PRINTER IS 16 !CRf 

3820 PRINT * * 'CLEAR CRT SCREEN 

3830 PRINT *CHECK TEHP. AT HOT AND COLO EDGE OF ADIABATIC ZONE tX-*-l/2La)* 
38X0 PRINT *1. IF TEHP. ARE O.X. TYPE i CONT TXPLOT* 

3850 PRINT * PRtSSt EXECUTE XEY* 

3860 PRINT *2. IF TEHP. NEED CORRECTING TYPEi CONT AOTEH * 

3870 PRINT ' PRESSi EXECUTE KEY* 

3880 PAUSE 'SEE ABOVE INSTRUCTIONS TO CONTINUE PROGRAH 

3890 RETURN 

3900 

3910 ' ALLOMS ONE TO CORRECT INPUT TLA AND TSA TEHPERATURES 

3920 AdtcMi INPUT ’TEHP. OF SAHPLE AT HOT END OF ADIABATIC ZONE(C)*,Tla 
3930 INPUT *TEnP. OF SAHPLE AT COLO END OF ADIABATIC ZOHECO’.Taa 
39X0 GOTO 360 


39J0 ' TXPLOT 

3970 ' PLOTS X vs T FOR -Li* < X < LiM 

3980 ' 

3990 Tiplot. PLOTTER IS 7,5,'9872A* 

XOOO DIH Titlet»(3) 

XOlO Tit letO(l)- Id* ‘SET PLOT ID. EQUAL TO PROGRAH 10. 

X020 ' INPUT ’MAIN TITLE TOR X vs T PLOT --18 CHARACTERS* , ll Uet* 1 1 ) 

X030 Titlat*<2)-*X(CM) 

XOXO Titlet»(3)-*TEHPTC>* 

X050 LIHIT 0,190,0,190 

X060 CALL AxeC-LlM,LlM.0,Th,l .S0,-Llm.0,l> 'IF LIH IS NOT AN INTEGER, AXE MILL A 
SSUHE THE NEXT LARGER IHTERGER 

X070 CALL LabIer(-LlM,Li«,0,Th,-LlM,0, TitlcttCD) 

X080 MOVE -LiM,Th 
X090 FOR I-l TO 51 

XlOO DRAW Xarray< I ) ,Tarray Cl ) 'DRAW X vs T ‘PLOT 51 PT. TEHP. PROFILE 
XllO NEXT I 

X120 ' DRAM Tmelt LlHE 

X130 HOVE -Liw.Tm 
XIXO DRAM LtM.TM 

X150 • DRAM LEFT SIDE OF ADIABATIC ZONE 

X160 Left--lat/2 
X170 HOVE Left ,0 
X180 ORAM Lef t ,Th 

X190 • DRAM HIDOLE OF ADIABATIC ZONE 

X200 MOVE 0,0 
X210 DRAM 0,Th 

X220 ' DRAM RIGHT SIDE OF ADIABATIC ZONE 

X230 MOVE -Lef 1 ,0 
X2X0 ORAM -Left . Th 

X250 GOTO Begin 'START ALL OVER 
X260 • 


•START ALL OVER 
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4?70 ' 

4280 • A4E 

4??S • SCALES AND DRAMS AXES 

4300 * 

4310 SUB AxeCXmtn.XiMi ,Ynin,V)Mll , Mint , Vint ,Xor 9 .Yorg , lypeJ ■ A>e 

Rev . A 


4320 DEE riSLow(Yor9-VMln\*.67»Het9htl 

4330 DEE EHLef t • (Xor 9 ~ Xmi n^> . 67*Mtdthi 

4340 CSIZE 2.5, .5 
4350 Helght'ViMi-Yiiiln 
4360 Mldth*XMax-Xiiiin 

4370 SCALE Xaln-ABSC C . 001 * . liFNLef 1 7*Midth} . X mi*ABS( C . 05* . lANOT FNLef tlMMldth 
) ,Y*ln-ABS<C. I*. 12«FHLow)NHelghl >,¥■•** ABSt(. 1*.12«M0T FHLow>BHet9ht> 

4380 CLIP Xnln , Xmax , Viein , Vaax 
4390 AXES Xml, Vint, Xorg.Vorg 

4400 DEC 
4410 LOIR 0 

4420 LORG 2 

4430 IF FNLeFt THEN LORG 8 

4440 Power-lNT<LOT(Vmt)) 

4450 LogacI-!NT((Type-l)/2) 

4460 Signal 
4470 Yend«Viiiax 
4480 FOR I-l TO 2 

4490 FOR Yy-Vorg TO Vend STEP SignRABS lY int > 

4500 IF CSign--l) AND lYyYorg) THEN Nxty 

4510 MOVE Xorg,¥y 

4520 IF (VyYorgi AND iXainOXorg) AND fXrwxv > X org) THEN MOVE Xor 


Nx 


g ,Vy*(F 
4530 
4540 
4550 
4560 
4570 
4580 
4590 
4600 
4610 
4620 
4630 
464 0 
4650 
4660 
4670 
4680 
4690 
4700 
4710 
4720 
4 730 
(FNLef t 
4740 
4750 
4 760 
4770 
4780 
4 790 
4800 
4810 
4820 
4830 
4840 
TMI H 


HLow-NOT FNLo«)«. OZaHeight 
Lab*Yy 

IF Log^cl THEN Lob-DROUNOl 1 0“Yy , 37 
GOSUB Label 

IF FNLeft then LABEL USING *# , * *Fal**' , X * i Lab 
IF NOT FNLeft THEN LABEL USING *# . X , * «Fat t iLab 

ty. NEXT Yy 
S ign- - 1 
Yend»Yain 
NEXT I 
LDIR 90 
LORG 2 

IF FNLow THEN LORG 8 
Power- INT<LGT<Xmt > ) 

Log4cl-NOT (Type MOO 2) 

S 1 gn- 1 
Xend-Xmax 
FOR 1-1 TO 2 

FOR Xx-Xorg TO Xend STEP SignaABSCXInt ) 

IF (Sign-l) AND <Xx-Xorg) THEN Nxtx 
MOVE Xx.Yorg 

IF CXx-Xorg> AND lYalnOVorg) AND (YaaxOYorgl THEN PLOT Xx* 
NOT FHLef t )•. 02»Width , Yorg 
Lab- Xx 

IF Logxcl Then Lab-DROUNO( lO^xx ,3) 

00SU8 Label 

IF FNLow THEN LABEL USING * « , * tFat •** . X * |Lab 
IF NOT FNLow IHEN LABEL USING * # , X . * «F at $ i Lab 

NEXT Xx 
Sign- - 1 
Xe.id- Xai n 
NEXT I 
OOT 0 Done 

bell IF (ABS(Lab)>-t00000) OR ( ABS (Lab) < . 001 ) OR (Pawer<--6) AND iLabOO) 
alt- 'MO. DE • 


Nx 1 1 1 
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4850 IF L«b*0 THEN Fntf'NO* 

4880 IF CABSCLab)>-100000) ON (ABSU.ab)< . 001 ) OR (P«»v*r<«-61 OR (Lab>0) THE 

H Ret 

4870 0l9>INTCLGT(ABSa«b})) 

4880 IF Logscl AND ( 0 l 9 < 0 ) THEN Fiat t- 'M. * OVALOtABStDlg) *1) 0*0* 

4880 IF 0ig>-0 THEN FntO-*H'«VAL»tDk9*l>0*0' 

4900 IF Logscl THEN Ret 

4910 IF (0ig<.0) AND (Dt 9 >Powcr> THEN FntO*‘n. *0VAL0<ABS(Pow«r)*l)0'0* 

4920 IF ( 0 t 9 < 0 } AND CDi 9 <-Power) 1HEN Fnt»-*M. ■•VALOtCABSiDlgl^DACOl 

9>-6)*6»(Dl9<>-6)>t*0* 

4930 IF (Po«er<0) AND (ABS(Po«ier X 6 -D 19 } AND «.Dlg>«0) THEN FatO-FattO 

*. *tVAL«CABS(Po«>er}*l>t*D* 

4940 IF (Power <0> AND (ABSCPower > >« 6 -Olgl AND (0lg>>0) THEN FwtO-FatO 

•*. •*VAL*(6-Dl9>*’0* 

4950 Kelt RETURN 
4960 Done I SUBEND 

4970 SUB Labler (Xatn , Kmk , Vain, Yatex , Xorg ,Vorg , Tit iet(*} ) ' Labler 

Rev. A 

4980 DEF FHLow*(Yorg-Valn(* .67aHeighli 

4990 DEF FNLef t-(Xor9-Xain<«.67a«mth) 

5000 Helght*ABS(YMx-yMlnl 
5010 Wldth>ABS(Xaax-Xmin) 

5020 LOIR 0 

5030 LORC 6 

5040 IF FHLow THEN LORO 4 

5050 HOVE Xaln»Wldtb/2,FNLowa(Vaax*.07aHeight}*NOT FHLo«»B(Yain-. 07aHelghU 
5060 LABEL USING */ ,K * 1 T 1 1 le» ( 1 J 
5070 LORG 4 

5080 IF FNLow THEN LORG 6 

5090 HOVE FNLef taXaax*NOT FNLef taxain. YorgrFNLowai . 05aHeight>-NOT FNLowRC.OSbH 
eight) 

5100 LABEL USING ' # , X • ( T 1 tie* 1 2) 

5110 LORG 4 

5120 IF FNLow THEN LORG 6 

5130 HOVE Xorg,FHLowa(YiMBx*.05aHei9nt)*N0T FNLowai Yai n- . OSBHe ight ) 

5140 LABEL USING *# ,X * , T i t let (3) 
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ABSTRACT 


The analysis of heat flow in directional solidification and < rystal 
growth by Bridgman- type processes has received considerable attention in 
recent years. Since the problon is complex, particularly when details! 
descriptions of furnace and sample properties are Included, numerical 
tech:iiques are often employed. While such techniques are Indispensable 
for modeling of real systems, it Is valuable to have available analytical 
methods, albeit simplified, to check and to guide the numerical analysis, 
to perform sensitivity analyses, and to perform trade studies on furnace 
design. 

This work develops a one-dimensional analvtlcal description of the 
heat flow in a translating rod that is applicable for Blot numbers less 
than unity. The model can accommodate different ?ropertles of the sam- 
ple in the solid and molten state, different heat transfer coefficients 
in the hot and cold zone, finite length effects In one zone, and approxi- 
mates ampoule effects by adjusting the Blot numbers to account for heat 
transfer in the ampoule. Also, the model contains an adiabatic zone M\d 
a booster heater zone that can be used to increase the thermal gradient 
near the solidification Interface or to provide an Independent control 
over the position of the Interface. Trade studies are performed on 
optimizing the length of the booster heater zone to obtain the maximum 
gradient In the sample without exceeding a specified maximum sample 


temperature . 


AM ANALYTICAL AfPROACH TO THERMAL 
MODELING OF BRIDGMAN-TYPE CRYSTAL GROWTH 

I. ONE-DIMENSIONAL ANALYSIS 

INTRODUCTION 

In Che design of vertlcsl Bridgaian-cype dlrectlonel solldlf Icetlon 
experiments, It Is Importsnt to be sble to predict the thermal profiles 
In the specliMn. Of particular Interest to the control of the process 
are the position and shape of the solidification Interface and the axial 
thermal gradients on each side of the interface. Generally, this problem 
requires the use of numerical methods to obtain the required solutions. 
While such techniques are necessary to take Into account the details of 
real systems, the Insight and ability to determine how variations In the 
controllable parameters affect the conditions near the Interface are 
diminished. For this purpose. It Is desirable to have an approximate 
analytical solution to guide the numerical modeling and the experimental 
procedures. 

Chang and Wilcox [1] solved the two-dimensional heat flow equation 
for a translating Infinite cylinder with uniform physical properties In 
a two-zone furnace. Equal heat transfer coefficients In the hot and cold 
zone were assumed in order to obtain an analytical solution. A one- 
dimensional analysis was used to investigate the shift in melt Interface 
position caused by the release of the latent heat of fusion and by finite 
length effects. These effects were subsequently Investigated by Rlquet 
and Durand [2], who obtained a restricted solution, and by Sukanek and Fu 
[3], who recently obtained a more general solution. Bartholomew and 
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HcXI»wtll uittd A c>n«*‘4Xmuifijlon«I Mtbcd t.o cbtJKin the 

taaipereturffi profile® in a dizftctiamlly »oli4ifl*4 A1 iSsMupLe, as a func- 
tion of Its position in the fumece, asBUTAing steedy~yt«t«{ conditions < 
Clyne [5] perforstcd both es:?ferlss«»tAl and coasputer analyaio of the 
directional solidiilcation of Al at higher trasimlatlonal vtlbcitiee and 
deiBOOstrated that for »et«rials with high thermal cninductivity there taay 
be a significant difference between the saaipl* pull rate and the actual 
grewth velocity. 

Fu and Wlico?! [6] invest i^ited a threc-roruj furnace in \d'dcr, the 
hot and cold rones were separated by a short adiabatic zone. Computer 
modeliag using a finite difference schem was employed, which allowed 
unequal heat-traasf «r coefficients to be. consildered in the hot and cold 
zones, It w»* eho*'^ that the use of m adiabatic zone significantly 
reduced the curvature of the leotherms between the hot and cold zones 
and allowed aruch better control over the. of the melt interface. 

This benefit, however, is at the exp«r..se of icrae loss in gradient. 

In many experiaoents involving saal.! sample diameters and fairly 
conductive saaterial;#, a oae-difflensicmal analysis is quite adequate for 
descr.ibing the conditions at the oelt interface, particularly if an 
adiabatic zone is used. Although such an analysis obviemsly cannot pro- 
vide InforciAticn about, the planarity of the interface, it can yield 
accurate predictions of the position of the interface and the jixial 
gradients in the solidification region as functions of tltc furnace and 
sample parameters. Thin provides a useful model for furnace optlalzation 
and for selection of experimental conditions. 
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The onc'^dimensionAl analyelt of Chang and Wilcox is sxtandad in 
this work to Include an adiabatic zone, different heat transfer coeffi- 
cients in the hot and cold zones, change in sample material properties 
associated with phase change, and the inclusion of a oooster iMater 
adjacent to the adiabatic zone to provide additional control over the 
gradient and position of the melt interface. 


ANALYSIS 


The heat-flow equation for a long, thin object with uniform cross- 
sectional area A and perimeter p and which moves along its axis with 
velocity u and transfers heat laterally with an external source or sink 
at temperature is [7] 



uoc 

k 


3X " Ak 


(T - T^) 


pc 3T 
k 3t 


0 , 


( 1 ) 


where K is the heat-transf er coefficient, k/ pc is the thermal diffusivity, 
and k is the thermal conductivity. For a cylindrical rod, the heat-trans- 
fer term becomes 


Ak a k 
o 


( 2 ) 


where a^ is the sample radius. 

It is convenient to introduce dimensionless coordinates, x ■ X/a . 

0 

The hear, transfer can be characterized by a dimensionless Biot number, 
Ha 


B. 

1 


k 


(3) 


and the motion of the sample can be described by a dimensiouless Peclet 
number 


P " ..ca u 
e 0 


k 


(4) 
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The dimensionless hest-flov equation for steady state beeches 

0 , ( 5 ) 




and the general solution Is 

T « T, -f C exp (a* x) + D exp (-a* x) 

where 


( 6 ) 


a* 


/ 2 2 
^UOL + P ± P 
e e 


and 


> 0 . 


If - 0, 

T - A + C exp (P x) 


( 7 ) 


To isolate various aspects of the problem, several simplified cases 
will be considered. First, the effect of the motion of the rod will 
be investigated in a three-zone furnace for a rod sufficiently long that 
end effects can be neglected. Next, the end effects will be Investigated 
for a motionless rod. Finally, the addition of a fourth zone, an 
Independently controlled booster heater between the main heater and the 
adiabatic zone, will be evaluated for its ability to increase the gradient 
in the sample at the melt interface and to control the position of the 
interface. 

Case I: Effect of Sample Translation 

The general solutions for an infinite rod in a three-zone furnace are: 


Tj^(x) - Tjj - D exp [S (x + Xj^)] 


. _ 00 < X < - X 


1 
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Tj^(x) ■ A + B exp 


Tg(x) - A* + B*exp (Ps x) 


- X, < X < X 

1 — — o 


X < X < X, 

o — — 1 


( 8 ) 


TgCx) - + C exp [-a (x - x^)] 


Xi ^ * 5. “ 


where 


/ 2 2 
/ + p/ - P, 


a* 


a 


2B, 




B* - 


P, + P- 
L L 


8 - rn 


H 


(9) 


subscripts L and S denote liquid and solid regions, respectively, and 

B„ and B^ are the Biot numbers in the hot and cold zones. Unless other- 
n C 

wise noted in the following, all lengths will be expressed in terms of 

sample radii. The length of the adiabatic tone, * 2Xj^. The position 

of the melt interface is x^ and is presuned to be in the adiabatic zone 

since this tends to yield a planar interface. 

The eqiiations must be solved simultaneously for seven unknowns. A, 

B, A*, B*, C, D and x^. Four constraints are provided by requiring T^ 

and T' to be continuous at -x.and T and T' to be continuous at x, . 

L 1 s s 1 

The remaining three conditions are provided by requiring T_ (x ) ■ T_(x ) 

L o So 


Tj^ and 




.-*x 


3Tj 

3x 


( 10 ) 


where is the latent heat of fusion term 


a u IH, 
L o f 


fi. H ""f 


(11) 


and is the heat of fusion. 
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and is the melting temperature. 

The position of the interface requires the solution of the tran- 
scendental equation 


\ ’'l ”h ■ “P 


exp ( 




+ L, 




I I + exp (Pj_x^) 




“P IJ - «p CPjX^) 


The remaining terms may be evaluated 


( 12 ) 


T - T 
H M 


exp (-P^x^)( — - lUexp (P^x^) 


( 13 ) 


B* 




“'P 1 1- exp (PgX^) 


( 14 ) 


A » Tjj + B exp (~Pl*i) 


>■) 


( 15 ) 


A* * T^ - S* 3xp (Pgxp I ^ + 1 


C = - 


PgS* exp (PgX^) 


( 16 ) 

( 17 ) 


P^B exp (-P^x^^ 


( 18 ) 


If Pg and P^ are small so that PgX^ ^ gradients in 


the sample at x mav be found from 
o 
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(T 




P, B e>:p (P^ X ) 

Li L 0 


H - ^ Vo ^ 


2 2 

* / * 1 ’'o ’'n 2 

x^ 1/B + - \[— - — 0(V 


(19) 


and 


G (x^) - P_ P* exp (P X ) 
so b So 


"m - V + '“s^o 


Xi^- 1/a - + P I ^ 


( 20 ) 

In the limit P « P ■ 0, the position of the melt interface can be 

1 j w 


obtained by equating kgG„(x^) - which yields 


X “ — 

o 


\ (Tr - ^)(^1 + 1/a) “ 

“h. «H- * '^S "c) 


( 21 ) 


and 


~ ’^M ~ 

a “ " + 1/a + 1/B 

A 


( 22 ) 


Case II: Investigation of End Effects 

To assess the importance of end effects, the complications intro- 
duced by the motion of the rod will be suppressed by setting the Peclet 
numbers to zero, and the rod will be inserted a finite distance L into 
the hot zone (see Figure 1). 
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The solution to the heat flow equation in the hot zone Is 


Tj^(x) - Tg - D exp [3(x+ x^)] + B exp [-9(x - x^)) ; -*3 1 x ^ -x^ . 

(23) 

In the adiabatic zone for P -^0, the heat-flow equation reduces to 
I'* - 0 and the solution is 


(x) ■ T + G- X ; -X, < X < X 

L o L 1 — — o 

and 


(24) 


* 

T (x) • T + G-X ; X ' X < X . 
o 0 b o 1 (25) 

As before, it is tacitly assumed that the interface is located in the 
adiabatic zone. 

The fourth equation is as in Case I, 


Ts(x) - T^ + C exp l-a(x - x^)] ; x^ < x < « 


(26) 


The interfacial conditions at -x^, x^ and x^ together with the 

requirement that T(x ) ■ T„ provide seven of the eight equations needed 

o n 

to specify the unknown coefficients. The final boundary condition Is 
imposed by the heat transfer to the end of the rod at -x^. 

--B* (Ty - T^(-X 3 )) . (27) 

H it it 

where B • H a /k and H is the heat-transf er coefficient between the 
rl o 

end of the rod and the furnace. 
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1 + B ^*/0 
1 - 


exp (-6L) 


exp [ti(L + L^)] -exp [-6(L - L^) ] 
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Case 111. Effect of Booster Heater 

The effects of an additional independently heated zone to sharpen 
the gradient In the sample and to provide additional control of the 
interface position will be investigated without the complications of 
sample motion or end effects. The heat flow eqiiations are; 


T (x ) - T jj - D exp 6 (x + x^) ; - * < x < -x^ 

T^ (x) » Tg + A exp [6(x + X2>] + B exp [-B(x + X2) ] ; -x^ ^ 

Tj_(x) -T^^ + G^x i -Xj<x<x^ (37) 

T (x) - T » + G X 
S o s 

Tg(x) - T^ + C exp [-a(x - x^)] 

where the length of the adiabatic zone Is ■■ 2x^, as before, and the 

booster heater is centered at -x and has length L (see Figure 2). 

/ B 

This set of equations contains nine unknowns. Interfacial boundary 
conditions at *1 supply eight of the necessary con- 
straining relations, and the requirement that supplies the 

ninth equation. Again, it is assumed that the Interface is located in 
the adiabatic zone. The solution of this set yields the following rela- 
tions : 


; X < X < X, 

o — — 1 

; X < » , 
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( T^ - I„) 11 -.xp (-34)1 -h (T„ - I^) t k^/k^_(Ty^ - I^) 
L. + 1 /B + 1/0 

A 

(k^/kj_) (I„-Ic> + C^ CXj + l/a> 


( 38 ) 


(39) 


‘'l 


(40) 


C .-Gj/O. 


(Al) 


B 


2 


exp (-BLg/2) 


(42) 


B exp (-ap/2) + G. /6 

A - 1 ® (43) 

exp (BLg/2) 

D - B exp (3L /2) - A exp (-6L /2) (44) 

O a 

T - T„ - G, x 
0 M L 0 

T * - T - x (45) 

o M S o 

Estimation of Effective Heat Transfer Coefficients 

To use the analytical model as a predictive tool, it is necessary 
to estimate the effective heat transfer coefficients in the hot and cold 
zones. Also, since the simplified analysis considers a homogeneous sample, 
some provision must be made to account for the presence of the ampoule. If 
the axial heat flow in the ampoule is very small compared with that conducted 
by the sample, the presence of the ampoule can be Included simply as an added 


resistance in the heat transfer between the sample and the furnace wall. A 


more general technique is to modify the thermal properties of the sample by 
defining an effective conductivity as the area weighted average of the sample 
and ampoule thermal conductivities. The effective Peclet numbers may be 
found by dividing the area weighted average heat carried by the sample plus 
ampoule by the effective conductivity. Since the ampoule does not contribute 
to the heat of fusion, the effective latent heat is the sample latent heat 
times the ratio of sample area to total area. 

At high temperatures, the heat transfer in the hot zone is dominated by 
radiation. The radiation heat transfer coefficient is given by 


H (x) - - 

^ad T^ - T^ix) 




(A6) 


where Q(x) is the heat flux transmitted at position x, is 


^e “ l/£* + a /a (1/e - 1) 

of f 


(47) 


e* is the effective emissivity of the ampoule /sample combination, e^ is 
the emissivity of the furnace wall, and a^ is the inside radius of the 
furnace muffle [8]. 

The heat transfer coefficient given by equation (46) is a function 

of position. The simplified one-dimensional analysis requires a constant 

value of that represents some type of average. The method chosen 

for computing this average is to require <H ,> to give the same total inte- 

rad 

grated heat transfer as H , (x) , i.e. , 
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Q(x) dx 



- Tj^(x)] dx 



- Tj^(x)^l dx. 


(A8) 


Since 


T^(x) - Ty - D exp [6* (x + x^)l. 


<H ,>• O F 
lad e 


T2 TI ■*■ IT \ ■*“ 4 \ 


(49) 


where T (-x ) ■ T„ - D and« is the temperature at the bo\mdary between 
LI H 

the hot zone and the adiabatic zone. This is not knoim ab initio but may 
be determined by an iterative process.^ 


In many cases the sample is enclosed in a quartz ampoule which has 
near unity transmission to some cutoff wavelength, X^, and is nearly opaque 
for X > X^. The effective emissivity of the ampoule /sample combination may 
be estimated using the technique suggested by Holland [9], 


e* - 0.95 [I - f(X )(l - E )1 (50) 

C 8 

where the 0.95 accounts for the surface reflectivity, E^ is the emissivity 

of the sample material inside the quartz ampoule for wavelengths shorter 

than X , and f(X ) is the fraction of radiation emitted at T„ which is 
c c n 

shorter than X . For metal-like samples, c is usually small and the long 
c s 

wavelength absorption in the quartz ampoule is the dominant mode of heat 
transfer into the sample. For such cases axial radiant heat transfer in 
the quartz ampoule is negligible. 


As a 


first estimate, (-x^ ) ^ 
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A siAllar traataent may ba employad In tha cold and by substituting 

T for and T (x, ) for T_(-x,). In this casa f(X ) is computad for 
c n s 1 L< i c 

s'*r‘ 

If thara is a fluid such as a lou-prassure gas in tha furnaca. its 
conductlva haat transfar is glvan by 


^cond a In (a,/a ) ’ 

O 10 

uhera k is tha conductivity of tha fluj.d. 

If flowing gas is usad for haat extraction, the haat transfer coef- 
ficient for forced convection is estimated from the work of Chen, Hawkins, 
and Solberg [10], 


H 

conv 


(X) 


, / \ 0.45 , ^ , 0 . 33 , , 0 , 8 , , 

(¥) (y (J) 


0.4 


(52) 


where v is the velocity, p is the density, is the heat capacity, y is 
the dynamic viscosity, and k is the conductivity of the gas; d is the 
equivalent diameter given by 2 (a^ - a^) , and x is tha position along 
the flow. Again, this heat transfer coefficient depends on position; and 
an average value must be found for use in the simplified model. This 
average is 


<H > 
conv 



H 


conv 


(x) 


dx 


H 

conv 


(L) 


.6 


(53) 


where H (L) is equation (52) evaluated at x - L, where L is the length 
conv 

of the sample in the furnace zone. 
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In the absence of forced convection, the heat transfer froio natural 
convection in a vertical furnace may be estimated by replacing v in equation 
(52) by v [g6 AT L]^ , where g is the acceleration of gravity, 6 is the 
coefficient of thermal expansion, and is the average temperature dif- 
ference between the ampoule and furnace wall. 

RESULTS 

The one-dimensional model described in this analysis was applied to a 
case that had been computed by a sophisticated two-dimensional heat- transfer 
analysis of a development furnace designed for use on the Space Shuttle as 
^<art of the NASA Materials Processing in Space program. The two-dimensional 
analysis was done with the System Improved Numerical Difference Analyzer code 
(SINDA) [10] which represents the 'umace by 580 nodes and 1420 conductors, 
and represents the sample plus ampoule by 5 nodes in the radial direction 
and 40 nodes in the axial direction. View factors are computed for radi- 
ation exchange between each element on the surface of the sample and all 
elements in the furnace within the field of view. 

The furnace consists of a 27.8 cm hot zc.:e, an indepeudently con- 
trolled 1.96 cm booster heater zone, a 2.03 cn adiabatic zone, and a 24.5 
cm cold rone. Heat extraction and thermal control in the cold zone are 
accomplished by flowing temperature-controlled He gas along the sample. 

The effective heat transfer coefficients and Blot numbers were esti- 
mated by the method outlined in the previous section. See the appendix 
for the numerical calculations. The thermal profiles of an infinite 

Mn-Bi/Bi eutectic sample were computed for T„ ■ 450*C, T_ ■ 40*C, and 

n C 

500'C. These are shown in Figure 3 together with similar 
temperature profiles computed using the SINDA code. Reasonably good 


agreement is obtained in the hot and cold zones » and the gradients at the 
solidification interface are in fairly close agreement. The major dis- 
crepancy between the two computation methods is in the position of the 
solidification interface, ihe interface locations predicted by SINDA are 
shifted approximately 3 tm closer to the cold zone than are those predicted 
by the one-dimensional model. There are several possible reasons for this. 


First* the simplified one-dimcnsitmal model assumes a perfectly 


adiabatic zone. The SINDA allows some heat transfer In the adiabatic zone 


from radiation losses and conduction from the hot (or booster heater) zone. 


Another possible reason for this shift is the fact that the simple model 
uses a constant average heat transfer cocfiicient in the hoc/booster 
heater zone. The inside diameter of the booster heater in the furnace 


in question is smaller than the diameter of the hot zone; hence, the 
conductive heat transfer will be underestimated in this zone by the 


simplified model. 

Despite these shortcomings, the simplified model can be used to 
investigate effects of varying sample and furnace parameters. For 
example, the effect of sample insertion length into the hot zone may be 
seen in Figure 4. As the sample is withdrawn from the hot zone, the 
interface shifts toward the heater. This effectively increases the 
solidification velocity above the pull rate. It may also be seen that 
the heat transfer coefficient into the «a6 of the sample makes only a 
small difference in the thermal profile and that the thermal field in 
the sample is almost fully developed after insertion lengths of a V 2 . * 

O Q 


"For finite samples whose "heat transfer length' 




as defined by Sukanek 


(v^ in Riquet and Durand notation) is » 1; this is equivalent to requiring 


the fraction immersed in the hot zone to be > 1/ 
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The effect of saaple translation it shown in Fisurs 5. As velocity 
is increased, the interface is shifted toward the cold zone. Note that 
the gradient in the Belt is less affected by translation chan the gradi- 
ent in the solid. Also note chat for the case of Mn-Bi/Bi eutectic, the 
release of latent heat of fusion is responsible for approxl&ately half 
of the interface shift. 

OptiuBUB Length for Booster Heater Zone 

It can be seen froa equation (38) chat the aagnitude of the gradient 

can be increased for a given by increasing However, if is aade 

long, it has virtually the saec effect as siaply increasing T . In aany 

H 

cases, the upper limit of T is fixed by other considerations, such as 

H 

melting point of the aapoule or vapor pressure of the sample, and teapera- 
cures higher chan I in the sample cannot be tolerated. The function of 

n 

the booster heater is to reduce the roll-off in temperature near the 
adiabatic zone without exceeding anywhere in the sample. This require- 
nent places a limitation on the length Chat is a function of the maximtsB 
T that can be supplied by the booster heater. 

O 

The requirement that T not be exceeded anywhere in the saaple can be 

n 

enforced by keeping D ^ 0. This implies from equation (44) that 

B exp (a„) > A exp (-GL^/2) . (54) 

The highest T. that can be used for a given L_ without overheating the sample 

D 0 

is obtained by setting D - 0 and using equatiens (42) and (43), 

'*3 - * -=L ' 

Ecuating this to (38), 
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(T T ) - ~ 

^ B Vmax ea^ + l/a) sinh BLg +"” :osh Tglg -~1 ' 

Substituting this into equation (55), the maximua gradient that can be 
achieved for a given Lg with no restriction on Tg is 


IG, i 


L ‘max 


3(L. + 1/a) s.^nh 6L_ + cosh 8L_ - i 
A D D 


(57) 


For BLg » 1, 




L ‘max 


+ 1/a + 1/6 


(58) 


which is the gradient given by eq. (22) for the three-zone configuta- 
tion. This comes about because eq. (56) forces T T„ as 6L » 1 to 

D n D 

avoid overheating the sample. On the other hand, for BL « 1, 

D 




L 'max 


+ 1/a + Lg/2 


(59) 


and the maximum gradient can be obtained by Lg 0. This, however, 
requires Tg -► ». Since Tj^^^x) fixed by practical considerations, 
there is an optimum length for the booster heat that depends on 
and the sample to be processed. This L (optimum) is the value of L 
for which the maximum gradient limited by heater temperature, eq. (33) 
Intersects the maximum gradient limited by sample overheating considera- 
tions, eq. (56). This value is 


^B(opt) 


1/8 in +/Ci^ + ^ , 


(60) 
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where 


1 - 3(1^ +1/0, - V 


(61) 


and 


L + i/a - 1 /e 

C m —n. .. - 

2 + 1 /a + 1/e 

A 


(62) 


The previous optimization assumed that the cold-end temperature 
remained constant. As T_ is increased, the solidification front x is 

D O 

driven toward the cold zone. If the cold-end performance permits, the 
T may be lowered to maintain the position of the interface near the 
center of the adiabatic zone. This also increases the maximum achievable 
gradient. 

From eq. (39), the interface can be maintained at x « 0 by setting 




(63) 


Putting this Into eq. (38) yields 


(Tg - Tg) [1 - exp (-pLg)] + 


+ 1/6 


(64) 


Employing eq. (55) to prevent sample overheating provides the 
maximum value for Tg - T^. 


(T. - T„). 


T - T 
H M 


B H'max Sx, sinh 6L + cosh 6L - 
I iS jd 


1 


( 65 ) 


Inserting this into equstlon (55) yields 
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max 


B(Ih - V 


sinh 6L^ 

slnh, + cosh 8^^ * ^ 


( 66 ) 


If L- Is less then the optimum length, the perfonumce Is limited by 

O 

booster heater performance and la given by eq. (64). As before, 
the optimum length for Is given by the Intersection of eqs. (64) and 
(66), which is 


S(opt) " 6 ^°*e 


1 + 0 + A + 20 + 


(I + ^)^ 


(67) 


where 


0 


^^B(max) 


-T)/(T - T) 

H'' ' '' H o'' 


( 68 ) 


The required cold>end temperature to position the interface at 
X ■■ 0 is obtained from equation (63) and becomes 


The effect of optimizing the booster-heater length is shown in 
Figure 6. In this case, the sample is Pb .Sn ,Te. It is 

• O • 

assumed that the booster heater cannot tolerate more than I200*C and 
chat the booster heater cannot be operated above 1600"C. As may be 
seen in Figure 6, the maximum allowable booster-heater temperature 
and gradient Increase as the booster-heater length is shortened until 
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1 


the maximum booster operating temperature is reached. Beyond this 
point, further decrease in booster heater length results In lower 
gradients as the booster heater loses its effectiveness. As would be 
expected, lowering the cold-end temperature to compensate for the 
booster heater in maintaining a stationary solidification Interface 
will give a higher gradient than operating with fixed T^. This can be 
seen more clearly in Figure 7, which compares the axial thermal profiles 
in a sample with no booster heater and I set to position the interface 
at X ■ 0, with profiles for an optimized booster heater operating a 1600 *C 
and fixed cold-end temperature and also with profiles for an optimized 
booster heater and T adjusted to position the melt interface in the 
center of the adiabatic zone. 

The booster heater provides another degree of freedom in position- 
ing the solidification interface independently of the hot and cold 
zone temperatures, which may be fixed by other considerations. Dif- 
ferentiating eq. (39) with respect to T^ yields 

^ 0 , 

(L^ + 1/a + 1/b) 

Since the sign of the derivative is positive, an increase in Tg will 
shift the melt interface toward the cold zone. The booster heater 
temperature may be selected to locate the solidification isotherm at 
a position of optimum planarity, or this temperature may be programmed 
to maintain the position of the isotherm as the sample is withdrawn 
from the furnace and to compensate for T„ changes as a result of com- 

Pi 

positional effects during solidification. 
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SUMMARY 


A relatively simple one-dimensional thermal model of the Bridgman 
growth process has been developed which is applicable to the growth of 
small diameter (£ 1 cm) samples with conductivities similar to those of 
metallic alloys. Although the model contains some idealized assianp cions 
that limit its ability to match actual thermal profiles exactly, it is 
useful for predicting thermal gradients, estimating solidification inter- 
face position, and analyzing effects of sample translation and sample 
insertion length. 

The fact that analytical results are available is particularly help- 
ful in performing furnace design trade studies and for inverting measured 
data to obtain furnace characteristics. Also, sensitivity analyses of 
furnace performaii , in terms of furnace and sample parameters may be 
accomplished by analytical techniques rather than by numerical analysis. 
This was used to determine the optimuQ length of the booster heater to 
produce a maxlmim thermal gradient in the sample for a given set of design 
constraints. 
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APPENDIX 


In the case analyzed by SINDA, T„ - 450®C, 1 • 40*C and 

H C 

■ 271.5. The Mn Bi/Bi eutectic was approximated thermally by pure 

Bi since this is a very low volume fraction eutectic. The thermal 

conductivities were taken as - 0.124 watt/cm/K and k^ ■ 0.072 watt/ 

cm/K. The quartz ampoule was 1.2 cm O.D. and 1.0 cm I.D. Thermal 

conductivity of the quartz was assumed to be 0.020 watt/cm/K. Effective 

conductivities are the area weighted average, or kj^(eff) ■ 0.0922 watt/ 

cm/K and k (eff) ■ 0,0561 watt/cm/K. The effective Peclet numbers are 

2 

given by (pc)eff a u/k(eff). Using p * 10.05 gm/cm for Bi and 2.3 

3 

gm/cm for quartz, c ■ 0.1463 J/gm for Bi and 0.493 J/gm for quartz, 
P^(eff) ■ 8.90 u and Pg(eff) • 14.63 u. Since the ampoule does not 
contribute to the heat of fusion, AH (eff) ■ 50.16 (.5/. 6)^ » 34.83 
Joules/gm. 

Cutoff wavelength for the quartz was taken to be 3.7 mm. For 
T„ - 723K, f(X ) - 0.20. and for T - 313K, f(X ) - 0.0015. The Z 

rl C C C S 

was taken to be 0.05. The effective sample/ampoule emissivity is there- 

"k k 

fore e « 0.7695 and e « 0.9486. Taking the inside diameter of 

H S 

tantalum muffle to be 5.2 cm with an £, » 0.3, the F term is 0.5441 for 

f ’ e 

the hot zone and 0.6279 for the cold zone. The average radiation heat 
transfer coefficients become 

<H^rad> « 0.00i45 watt/cm^/K 

2 

<H^rad> - 0.000557 watt/cn K . 


Heat extraction In the cold zone la accomplished by flowing He. 
Using a 4 ib/hr flow rate through the cold zone (a^ * 1.77 cm, a^ - 

-0 4 

0.6 cm), the heat transfer coefficient used in SINDA is 23.36 x(ft) 

2 

BTU/hr/ft /*F. For a 25 cm long sample inserted half its length into 
the cold zone, the average heat transfer coefficient is 58.9 BTU/hr/ 
ft^/"F or 0.0312 watt/cm^/K. 

A stagnant region of He will also be present in the hot end. 

Since the furnace was designed for low-gravity operation, convective 

heat transfer was not analyzed. The conductivity of He at 450*C is 

taken as 0.0023 watts/cm/K, and the conductive heat transfer coefficient 

2 

is therefore 0.00261 watt/cm /K. 

The total effective heat transfer coefficient is the svm of H(rad) 
H(cond) and H(conv) . The effective Biot number is the total effective 
heat transfer coefficient multiplied by the ampoule radius divided by 
the effective conductivity. The values are: 

Hjj - 0.00706 watt/cm^/K, - 0.04594 

H ■ 0.03176 watt/cm^/K, B ■ 0.3396 

c c 
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Figure 2. 


Figure 3. 


Figure 4. 


Figure 5. 


Figure 6. 


Figure 7. 


Furnace /sample configuration shoving coordinate system used to 
investigate end effects. 

Furnace/sample configuration showing coordinate system used to 
investigate the effects of a booster heater. 

Comparison of thermal profiles for a stationary Mn Bi/Bi eutectic 

sample computed with the simplified one-dimensional model with 

those computed by the SINDA code for T_ ■ T„ and T_ • T„ + 50®C. 

on on 

Effect of sample insertion length on thermal profile for a 
stationary Mn Bi/Bi sample in a furnace with T - . 

D n 

Effect of sample translation of thermal profiles for an infinite 
sample moving at 0, 10, 20, 30 cm/hr in a furnace with . 

Effect of booster heater length on maximxan gradient obtainable in 
Pb pSn _Te sample without exceeding sample temperature 
of 1200®C. As heater length is shortened, the allowable booster 
heater temperature is Increased (with a concomitant increase in 
thermal gradient) until the upper limit imposed by heater materi- 
als is reached (assumed in this case to be 1600*C). At shorter 
lengths, the sample gradient is limited by booster heater perform- 
ance. The optimum length for the booster hsater is the intersection 
of these two curves. 

Thermal profiles in Pb j,Sn -Te with optimized booster 

heater analyzed in Figure 3. T^j^ was chosen to position the 

solidification interface at x *0 for T_ ■ T. 
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Figure Cap t lone 


Figure I. 


Figure 2. 


Figure 3 


figure 4. 


Figure 5. 


Figure 6. 


Figure 7. 


Fumace/aaffiple configuration showing coordinate system used to 
investigate end effects. 

Furnace /sample configuration showing coordinate system used to 
investigate the effects of a booster heater. 

Comparison of thermal profiles for a stationary Mn Bl/Bi eutectic 

sample computed with the simplified one*^imenslonal model with 

« 

chose computed by the SINDA code for T_ * T„ and T_ ■ T„ 50*C. 
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Effect of sample insertion length on thermal profile for a 
stationary Mn Bi/Bi sample in a furnace with T - T . 

0 n 

Fifecc of sample translation of thermal profiles for an Infinite 
sample moving at 0, 10, 20, 30 cm/hr in a furnace with . 

Effect of booster heater length on maximm gradient obtainable in 
Fb gSn >Te sample without exceeding sample temperature 
of 1200*C. As heater length is shortened, the allowable booster 
heater temperature is increased (with a concomitant increase in 
thermal gradient) until the upper limit Imposed by heater materi- 
als is reached (assxmed in this case to be 1600*C). At shorter 
lengths, the sample gradient is limited by booster heater perform- 
ance. The optimmn length for the booster heater is the intersection 
of these two curves. 

Theraal profiles in Pb „Srs »Te with optimized booster 
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heater analyzed In Figure 3. was chosen to position Che 

solidification interface at x - 0 for T_ • 7„. 
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